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DATA MANAGEMENT REQUIREMENTS STUDY 
FINAL REPORT 

INTRODUCTION 

The AMPS Data Management Requirements Study has been aimed at defining 
and developing the initial data simulation of AMPS Instruments and associated 
control and display functions required to perform controlled active experi- 
ments from the AMPS laboratory. The four tasks that were undertaken to ful- 
fill this aim are as follows: 

- Develop and provide to NASA the data output (raw and formated) 
for specific sensors and detectors. 

- Provide the signal and information formating required to activate 
and control AMPS Instruments and Instrument support units. 

“• Specify display data and display formats. Provide supporting 
flowcharts for the AMPS simulation 

- Develop requirements and procedures for simulating a representative 
set of AMPS experiments. 

This report has been organized to provide a comprehensive users guide for 
the data requirements and software developed for the following experiments: 

1. Electromagnetic Wave Transmission Experiment. 

2. Passive Observation of Ambient Plasmas. 

3. Ionospheric Measurements with a Subsatellite. 

A. Electron Accelerator Beam Measurements. 

5. Lidar Trace of Acoustic Gravity Waves in the Sodium Layer. 

In addition, the data requirements of a sixth experiment have been 
developed wi thou t the associated software flow charts. This experiment is 

6. Determination of the Wake of a Test Body. 

The data output of each of the sensors and detectors are presented in 
Section 1.0. This section essentially satisfies the requirements of Task 1 
described previously. Section 2.0 defines the control requirements called 
for in Task 2. The complete simulation of 6 experiments is given in Sec- 
tion 3.0 of this report. The software requirements and flow charts for five 
of the six experiments described in Section 3*0 are presented in the Appendix. 


Task 1 
Task 2 
Task 3 
Task A 
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1 .0 RAW AND FORMATED DATA OUTPUT FOR 
AMPS SENSORS AND DETECTORS 

The output requirements for the sensors and detectors listed in Table 

1.0 are described in Section 1.1. The formated data is given in Section 1.2. 

1.1 OUTPUT REQUIREMENTS 

1.1.1 Remote Subsatellite Systems (2) 

1.1. 1.1 TV System. The subsatellite television system will be used 
to observe optical emissions at large distances from the Shuttle, its data 
output will include its pointing direction, its sensitivity, the spectral 
content and finally the television picture itself. Because the data rate 
for the optical imaging will be so high, >10^ bps, the data must be handled 
by special video channels and will not be processed by the spacelab computer. 
The pointing data, the spectral content (if a filter system is used) and the 
sensi ti vi ty data are all moderately slow and require only one 16 bit word 
once per second. Hereafter, housekeeping data will not be speci f ical ly de- 
scribed. However, when the telemetry frame format is presented there will 

be enough room in the format to accomodate all of the subcommutated house- 
keeping data. 


Table 1.0 Sensors and Detectors 
A. On Two Remote Subsatellite Systems 


1. 

TV System 

6, 

Ion Mass Spectrometer 


2. 

Triaxlal Fluxgate Magnetometer 

7. 

Triaxial Hemispherical 

Analyzer 

3. 

Search Coil Magnetometer 

8. 

VLF/RF Receiver 


k. 

Cylindrical Electron Probe 

9. 

Electric Field Meter 


5. 

Segmented Planar Trap 




On 

Diagnostic Boom A 




1. 

One Meter Loop Antenna 

8. 

ion Mass Spectrometer 


2. 

Short Electric Dipole Antenna 

9. 

Spherical Ion Probe 


3. 

Tri axial Search Coil Mag- 

10. 

Cylindrical Electron Probe 


netometer 

11 . 

Planar Segmented Probe 


4. 

Rubidium Magnetometer 

12. 

Neutral Mass Spectrometer 

5. 

Triaxial Fluxgate Magnetometer 

13. 

Triaxial Hemispherical 

Analyzer 

6. 

33 Meter Electric Dipole Antenna 

14. 

Planar Electron Trap 


7. 

A1 ignment TV 




On 

the Remote Sensing Platform 




1 . 

Scanning Grating Spectrometer 

3. 

Filter Photometers 


2. 

Fabry-Perot Intferferometer 

4. 

Total Energy Detector 
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1.1. 1.2 Trl axial Fluxgate Magnetometer . This magnetometer will be 

be used to measure the magnetic attitude of the subsatellites and to measure 
slow period magnetic waves. The three axes will be mutually orthogonal. 

The maximum ambient signal could correspond to 50,000 gammas either positive 
or negative and a sensitivity of at least 1 gamma will be needed. Therefore 
there is a need for three channels of data each able to handle at least a 
range of 100,000 with accuracy of 1 part in 100,000. While the fluxgate 
data is actually analog it is assumed that the magnetometer will digitize 
its own data. The digital data format of the fluxgate requires 17 bits for 
each of the channels. A sampling rate of 200 times per second will be 
needed to obtain the proper number of points to do magnetic wave analysis. 
The fluxgate data will fit into four 16 bit words with 17 bits for each axis 
3 bits for polarity and 10 bits for housekeeping. 

1.1. 1.3 Search Coil Magnetometer . The search coil magnetometer will 
be used to detect magnetic field variations at frequencies above 30 Hz up 

to a limit Of perhaps 10 or 20 MHz. The magnetometer should be able to give 
wideband data for this whole range of frequencies which means that it will 
exceed the capabilities of the spacelab computer and must be handled 
like a video signal. The data consists of a voltage signal level with a 
20 megahertz bandwidth and a dynamic range of as much as 10^. 


1.1. 1.4. Cylindrical Electron Probe. Cylindrical electron probes 
come in many varieties. The one used for this simulation measures ambient 
electron densities and temperatures by applying a sweep voltage to the 
probe. The data will consist of an electrical current value and a vol tage 
value. The data are assumed to be digitized by the probe's own electronics. 
The current wi 1 1 take 20 bits and the voltage 10 bits each sampled at a rate 
of 1,000 times per second. The data format consists of two l6-bit vvords. 


1.1.1. 5. Segmented Planar Trap. The segmented trap is used to mea- 
sure electron or ion densities, temperatures, or dri f t veloci ti es . Usually 
it is used in the ion drift velocity mode of operation. The measurement 
consists in comparing the electrical current values for each of the segments 
For this simulation we will assume that the trap has four separate segments 
and that there will also be a retarding potential voltage. The dynamic 
range of the current can be quite large, therefore each of the four segments 
of the trap will be given 20 data bits and a time resolution of I/IOOO 


seconds. The retarding potential voltage wM 1 get II data bits at the same 
time resolution. These data require six l6-bit words 1,000 times per second. 

1.1.1. 6 Ion Mass Spectrometer. The Ion mass spectrometer will be 
used to separate the various ionospheric ions by mass and to determine either 
the relative or absolute ion density distributions. If only the ambient 
plasma Ions were to be considered, a dynamic range for the density of 10° 
would be needed, however during some AMPS missions artificial plasma clouds 
will be made and their densities could be a factor of ten larger than ambient 
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raising the needed dynamic range to 10 . The mass resolution of the spectrom 
eter needs to be at least one part in forty and the masses measured could 
possibly be as large as several hundred AMU for artificially produced Ions. 
The sampling rate will be AOO per second. The digital format consists of 
two l6-bit words, 21 bits for the ion density current, 9 bits for the mass 
and 2 polarity bits. 

1.1. 1.7 Triaxial Hemispherical Analyzer. The hemispherical analyzers 
will be used to measure the flux and energy of charged particles whose 
energy lies in the range from 50 to A0,000 electron volts. The energy spec- 
trum of the particles Is obtained by sweeping a voltage between the two 
analyzer plates, and both electrons and ions can be measured by changing the 
polarity of the sweeping voltage. The maximum counting rates will be on the 
order of 10^ particles per second. A fast time resolution will be needed 
because of the complex temporal behaivor of energetic particles. The tri- 
axial analyzer needs three separate detectors, but only one sweeping voltage 
which is applied to the three analyzers. Three 16-bit words with 13 bits 
devoted to each of the three counting data channels and 9 bits. for the sweep 
voltage. These three words wi 1 1 be sampled 5,000 times per second. 

1.1.1. 8 VLF/RF Receiver. The VLF/RF receiver will be used to measure 
the amplitude of electric field fluctuations in the frequency range from 300 
Hz to 20 MHz. The data will be contained within a 20 MHz broadband channel . 
The dynamic range may be as large as 10^. Because the bandwidth will be so 
large these data must be treated like the television and search coil mag- 
netometer data and processed outside of the space! ab computer. 

1.1.1. 9 Electric Field Meter. The electric field meter will be used 
to determine the amplitude and direction of the electric field vector in the 
frequency range from 0 to about 30 Hz. It will be a three axis detector. 
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Each axis has as an output a signal voltage level which is proportional to 
the electric potential difference between two sensors. The dynamic range 
of the instrument is about 10^ and the sampling rate should be about 300 per 
second. This requires three 16-bit words with 15 bits for the signal voltage 
and 3 bits for polarity, 

1.1.2 Diagnostic Boom A 

1.1. 2.1 One Meter Loop Antenna. The one meter loop antenna will be 
used to make measurements of higher frequency magnetic field variation. The 
frequency range will depend upon design details, but should cover the range 
from 1000 Hz to 20 MHz with a dynamic range as large as 10^. The data chan- 
nel must have a 20 MHz bandwidth and can not be processed through the space- 
lab computer. These data will be treated like video data. 

1.1. 2. 2 Short Electric Dipole Antenna. The short electric dipole 

antenna is sensitive to electric field variations at the same frequencies 

as the loop antenna, from 1000 Hz to 20 MHz. Its dynamic range can be as 

large as 10°. Its data can not be processed through the spacelab computer 
because of the large bandwidth, but must be treated like a video signal. 

1.1. 2. 3 Triaxial Search Coil Magnetometer. Each axis of this instru- 
ment will be exactly like the subsatellite search coil magnetometer treated 

earlier. With three axes the bandwidth will be increased to 60 MHz and the 

data must be processed on video channels. 

1.1. 2. 4 Rubidium Magnetometer. The rubidium magnetometer will make 
accurate measurements of the value of thei total magnetic field at its loca- 
tion. This Instrument is capable of accuracy of from 1 to 10 parts in 10^. 
This magnetometer needs 23 bits for signal level and one bit for polarity 
200 times per second to read out its data. This will require two l6-bit 
words and will leave 8 bits to be used for housekeeping. 

1.1. 2. 5 Triaxial Fluxgate Magnetometer. This instrument has been 
described previously. 

1.1. 2. 6 33 Meter Electric Dipole Antenna. The 33 meter dipole antenna 
measures VLF/RF electric fields in the frequency range from 300 Hz to 10 
megahertz. This antenna can also be used as a low power transmitting an- 
tenna in the same frequency range. When used as a receiver its data output 
is the same as the VLF/RLF receiver. 
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1.1 *2.7 Alignment Television. The alignment television system will 
be used to visually determine the position and orientation of the platform 
at the end of boom A and also to act as a viewing system for experiments 
carried out close to the Shuttle. If the system is used only for alignment, 
it can be black and white. But some experimental observations will need an 
accurate color system provided by an optical filter wheel. The data from 
the alignment television system will be assigned to a video channel and can 
not be processed through the spacelab computer. 

1.1. 2. 8 ion Mass Spectrometer. This ion mass spectrometer will be 
identical to the one described in Section 1.1. 1.6. 

1.1.2. 9 Spherical Ion Probe. The spherical ion probe will be used 

to determine the ion density and temperature and will give some data on the 

distribution of the ion species. This probe operates by sweeping the probe 

voltage and measuring the probe current. The dynamic range of the current 
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can be as large as 10 . This requires 20 bits for the current while 10 bits 
for the sweep voltage on 2 polarity bits will also be needed. Two 16-bit 
words read 1,000 times per second will be sufficient for the spherical ion 
probe. 

1.1.2.10 Cylindrical Electron Probe. This cylindrical electron probe 
will be identical to one on the subsatellite. 

1.1.2.11 Planar Segmented Probe. The planar segmented probe will be 
identical to the segmented planar trap described earlier. 

1.1.2.12 Triaxial Hemispherical Analyzer. This trlaxlal hemispherical 
analyzer will be identical to the one described earlier. 

1.1.2.13 Planar Electron Trap. The planar electron trap is used as 
a retarding potential analyzer. The trap measures the current of electrons 
striking it as a function of the swept retarding potential and allows the 
electron density distribution function to be calculated. The current measure 
ment will have a dynamic range of 10^ requiring 20 data bits. The voltage 
will require 11 data bits. Two 16-bit words with a repetition rate of 1,000 
times per second will be sufficient for the planar electron trap data. 

1.1.3 Remote Sensing Platform 

1.1. 3*1 Scanning Grating Spectrometer. The scanning grating spectrom- 
eter will be used to obtain moderate wavelength resolution of atmospheric 
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optical emissions over the range from 250 to 1,000 nanometers. The photo- 
detector will be used in both a counting and a current mode with a dynamic 
range of 10^ for the count mode and 10^ for the current mode. The wave- 
length resolution should be about 0.01 nanometers and would require a range 
5 . 

of 10 . These requirements will give 2M data bits for the counts, 17 data 
bits for the current and 17 data bits for the wavelength. The sampling 
rate will be 1,000 times per second. Four 16-bit words could be used with 
6 bits left over for housekeeping functions. 

1.1. 3*2 Fabry-Perot Interferometer. The Fabry-Perot interferometer 
will be used to obtain very high resolution spectra of atmospheric optical 
emissions in a wavelength region from 300 to 2,000 nanometers. The dynamic 
range for counting photons will be about 10^ and for wavelength resolution 
will be about 10^. The time resolution and the interferometer scanning 
time require 1,000 data samples per second. Three 16-bit words can be used 
to handle the Fabry-Perot interferometer data with 2^ data bits for wave- 
length, 17 data bits for counts and 7 data bits for housekeeping. 

1,1. 3*3 Filter Photometers. The AMPS Shuttle payload will carry as 
many as eight photometers with separate filters on each used to delineate 
optical lines in the atmospheric emissions. The photometers will be used 
in the photon counting mode. They will be sampled 1,000 times per second 
and wM 1 be a 1 located one 1 6-bit word for each photometer. 


1.1. 3*^ Total Energy Detector. The total energy detector will mea- 
sure the solar energy incident on the Earth in several wavelength intervals. 
This instrument does not need rapid time resolution nor does it need a wide 
dynamic range within a given Wavelength interval. By using subcommutation 
within the instrument all of the total energy detector data can be handled 
by one 16-bit word read 100 times per second. 

1.2 DATA FORMATING 


The data formats for the sensors are given in Figure 1.0 through 1.5. 
There are many al ternative ways of presenting the data, and the format se- 
lected is sufficiently general so that it can match the requirements of any 
data management system. The information flow rates shown in the figures 
range from 200 to 5000 measurements per second per instrument. Although the 
higher rates are beyond the capability of the CVT computer to generate and 
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process, they represent what would be desirable from an experimenter's view- 
point if Unlimited resources were available. It is noted that in the com- 
puter flowcharts for the CVT simulation presented in the Appendix the highest 
sensor data rates called for are about 100 measurements per second (lidar, 
Experiment #5) » which is well within the capabilities of the Spacelab com- 
puter and CVT simulation computers. 

1.2.1 Remote Subsatell i tes (2) 

In this simulatiori the two subsatellite data requirements are treated 
as being Identical. The discussion below encompasses the needs of one sub- 
satellite. 

The instruments on the subsatellite require both broadband analog and 
digital data transmission to the Shuttle. The television system, the search 
coil magnetometer and the VLF/RF receiver all together require a 60 mega- 
hertz broadband analog transmission which must have a very large dynamic 
range on a frequency modulated channel. 

The digital data has been assumed to be grouped into 16-bit words in 
order to be compatible with both the spacelab computer and the Sum C system. 
There is a requirement for 25,500 words per second from the subsatellite 
instruments. The characteristics of the words for each of the instruments 
Is shown in Figure 1.0. A possible telemetry format for these digital data 
uses a 16 by 16 word frame repeated at a rate of 100 times per second. Frame 

identifiers arid housekeepi ng data are allocated to the first word In the 

! 

frame and the other 255 words are comprised of the instrument words. Besides 
the first l6-bit word there are 16 bits per full frame available in the 
fluxgate magnetometer words for housekeeping data for a total of about 3,200 
bits per second. Proper subcommutation of these bits should provide ade- 
quate Instrumental housekeeping data. The housekeeping for the other sub- 
satellite systems will require their own telemetry allocations. Figure 1.1 
shows a schematic table of a possible instrumental main frame format for 
the subsatellite instruments. 

1.2.2 Diagnostic Boom A 

Because of the maneuverability Of the boom It will probably be neces- 
sary to either telemeter the instrument data to the spacelab or have at most 
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Figure 1.0. Possible Allocation of Digital Words for Subsatellite Instruments 
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Figure 1.1 Possible Telemetry Format for the Subsatellite Instruments 




11 


one or two very wide band coaxial cables. Therefore there is a need to for- 
mat the transmission of data from the diagnostic boom. 

The loop antenna, the short electric dipole antenna, the triaxial 
search coil magnetometer, the 33 meter dipole antenna, and the alignment TV 
all require wideband analog data channels. The total bandwidth for all of 
these instruments is \k0 megahertz and must be accomodated on frequency 
modulated signals. 

The digital data requirements of the diagnostic boom instruments are 
29,300 1 6-bit words per second. The instruments on the boom are identical 
with those on the subsatellite with the following exceptions: (1) there is 

no electric field meter on the boom, (2) a spherical ion probe is added to 
the boom, (3) a rubidium magnetometer is added to the boom, and (4) a 
planar electron trap is added to the boom. Figure 1.2 shows the word 

formats for these three added instruments. A possible telemetry format for 
these digital words uses a I6 by I9 word frame repeated 100 times per second. 
There are 13 housekeeping words, which may also be considered to be spare 
data channels for instruments that might be added to the diagnostic package, 
as well as 16 bits per frame in the fluxgate magnetometer and I6 bits per 
frame in the rubidium magnetometer words. Figure 1. 3 shows a possible con- 
figuration of this main telemetry frame for the diagnostic boom A. 

1.2.3 Remote Sensing Platform 

The data channels from the instruments on the remote sensing platform 
might be hardwired from each Instrument to the spacelab data bus, however 
the interface might be simpler if the data were mutliplexed and sent on one 
signal wire. In that case the data from the remote sensing platform require 
15»100 words per second. Figure 1.^ shows a possible format for data from 
each of the separate instruments. Figure 1.5 shows a possible format for 
the main data frame using a 16 by 10 word array with 9 words for housekeeping 
or as spares. There are also 60 bits in the scanning grating words and 70 
bits in the Fabry-Perot words for further housekeeping data. 


Rub i d i urn Magnetometer 
2 Words - 200 Times per 
Second 


Spherical Ion Probe 
2 Words - 1 000 T i mes per 
Second 


Planar Electron Trap 
2 Words “1000 Times per 
Second 



Figure 1 .2. Possible A1 location of Digital Words for Diagnostic Boom A Instruments 
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H * Housekeeping + Identifiers + Spares 
F = Fluxgate Magnetometer 
C = Cylindrical Probe 
S * Segmented Planar Probe 
I “ Ion Mass Spectrometer 
A * Hemispherical Analyzer 
R = Rubidium Magnetometer 
P = Spherical Ion Probe 

T = Planar Electron Trap. Frame is repeated 100 times per second. 
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Figure 1.3. Possible Frame Format for Diagnostic Boom A Instruments 
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Figure 1.4. Possi ble Allocation of Digital Words for Remote Sensing Platform Instruments 
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13 14 15 16 


G = Scanning Grating Spectrometer 
FP = FabrayPerot Interferometer 
PI = 1st Filter Photometer 
P2 - 2nd Filter Photometer 
P3 3rd Filter Photometer 
P4 = 4th Filter Photometer 
P5 = 5th Fi Iter Photometer 
P6 = 6th Fi Iter Photometer 
P7 = 7th Filter Photometer 
P8 = 8th Fi 1 ter Photometer 
TE = Total Energy 

Frame is read 100 times per 
second from left to right, 
top to bottom. 
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Figure 1.5. Possible Frame Format for Remote Sensing Platform Instruments 
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2.0 SIGNAL AND INFORMATION FORMATING REQUIRED 
TO -CONTROL THE AMPS INSTRUMENTS AND 
INSTRUMENT SUPPORT UNITS 

The control requirements for the instruments and instrument support 
units listed in Table 2.0 are described below. In each case the first 
approximation to the information formating is supplied in the form of a 
conceptual block diagram showing the interface between the instrument sys- 
tem, the control buffer of the spacelab computer and the input-output (dis- 
play) data* 


Table 2.0 Instrument and Instrument Support Units 


A. Remote Sensing Platform 



B. Accelerator System 



1. 

Ion Accelerator 

4. 

The Gimbaling System 

2. 

Electron Accelerator 

5. 

The Accelerator Pointing System 

3. 

MPD Arc 

6. 

Accelerator Deployment System 

C. High Power Transmitter/Antenna System 


1. 

10 kw, 200 kHz - 2 MHz 

3. 

1 kw, 300 Hz - 200 kHz 


Transmitter/Coupler 


Transmitter/Coupler 

2, 

10 kw, 2 MHz - 20 MHz 

4. 

1000 Foot Dipole Elements 


Transmitter/Coupler 

5. 

Transmitter Deployment System 

0. Diagnostic Boom A 



1. 

Boom A Deployment System 

5. 

Electric Dipole Extension 

2. 

Boom A Maneuvering Functions 

6. 

Power Supply 

3. 

Gimbaled Platform 

7. 

Data System' 

4. 

5 Meter Sub -booms 

8. 

Alignment TV System 

E. Boom B 



1. 

Boom B Deployment System 

4. 

Wave Generator 

2. 

Boom B Maneuvering Functions 

5. 

Low Energy Electron Gun 


3. Target Deployment 6 Retraction 

F. Deployable Units 

G. Deployable SubsatelHtes 

1. Deployment Mechanism k. Telemetry System 

2. TV System 5. Ranging and Control System 

3. Transponder 


17 


2.1 REMOTE SENSING PLATFORM 

The control requirements for the remote sensing platform derive from 
the fact that the platform's pointing direction must be controlled so that 
the optical instruments placed on the platform can be pointed toward spe- 
cific targets or regions of space, in some cases the platform will be 
turned to a Specific direction and be held stationary, in other cases the 
platform will execute a scanning mption either in response to a prepro- 
gramed set of Instructions or to an experimenter's active control commands. 

In all of these cases the Initial commands will be given by the experimenter 
and the spacelab computer will serve as a buffer between the experimenter 
and the platform control motors. Figure 2.0 shows a conceptual blocJk dia- 
gram of this command-con trol-d I splay scheme. 

The pointing accuracy requirements on the remote sensing platform 
will be on the order of arc min relative to a coordinate system fixed 
on the Earth's horizon. For this reason It is likely that the position 
sensing unit will probably be referenced against an Inertial coordinate 
generator. During either the point-and-hold mode or the preprogramed scan 
mode the computer will be' in full control of the platform and will display 
position data to the experimenter. During the experimenter scan mode, the 
computer will act as a buffer and will do position calculations for display. 

The control format for the experimenter or the programed scan consists 
of specifying two angular coordinates with an accuracy of better than 1 arc 
per minute. The experimenter will simply type in the two desired angles 
and the computer will translate the command into a motor control signal. 

The display output will be a realtime reading of the two angular values. 

For the purposes of a simulation of the remote sensing platform all one 
needs to specify is the rate at which the motors run and the time at which 
they begin running. The experimenter or the scan program will set the 
position to which the platform is to move and will start the platform motors 
The computer will then run the angle variables as a function of time until 
the preset values are reached, and will display these values continuously. 
These angle values will also be needed to simulate the output of the instru- 
ments mounted on the platform. In actual use on board the Shuttle the 
computer will compare the preset angle values with the platform position 
sensors and turn the motors on and off In a manner to bring actual posi- 
tion into agreement with preset position. 



Figure 2.0 Conceptual Block Diagram of Remote Sensing Platform Control 
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2.2 GIMBALED ACCELERATOR SYSTEM 

As the AMPS accelerator system is currently planned, there will be 
no gimbaled mount of the accelerator system. The directioning of the beams 
will be accomplished with magnetic bending coils. Therefore, in the fol ' 
lowing descriptions the bending coils rather than the gimbal system will be 
used. 

2.2.1 Ion Accelerator 

The controls for the ion accelerator involve the regulation of the 
ion source^ the choice of ions, the accelerating voltage, the beam current 
and the neutralization of the beam. Because the cyclotorn radius of even 
the hydrogen ion is large compared to the dimensions of the ion accelerator, 
no direct means is provided for varying the beam ejection direction. If the 
pitch angle of the ions must be varied, the Shuttle attitude must be changed. 
The control functions divide Into two categories: accelerator setup and 

beam characteristics. Figure 2.1 shows a conceptual block diagram of the 
ion accelerator control scheme. 

The setup of the ion accelerator requires that the source power supply 
be turned on, that the current and voltage to the ion source be adjusted, 
that the gas supply to the source be selected for ion type and the flow 
regulated, and that the neutralizer power supply be turned on and the 
neutralizer be heated hot enough for copious electron emission. The on-off 
commands will be simple level changes initiated from the keyboard. The 
regulation of the source current and voltage and the neutralizer current and 
voltage requires that the experimenter be able to see a realtime display of 
these quantities while he selects new values. For the purposes of simulation 
the voltage and current values will be identical with the selected values. 
During an actual flight, there wi 1 1 be variations between selected and actual 
values. The regulation of the gas flow rate requires a realtime display of 
the source pressure whi le the experimenter selects new values of the flow 
rate. In the simulation there wi 1 1 be required a functional relationship 
between the source pressure and the flow rate while during an actual flight 
the pressure wi 1 1 be displayed as an independent variable. The selection 
for type of gas, which determines the ion species of the beam, will be 
given by the experimenter on the keyboard. 
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The control of the beam characteristics will be done either by a pre- 
planned program or by the experimenter from the keyboard. Only the beam 
current and accelerating potential will be controlled. Both the accelerator 
voltage and the current are controlled by a voltage signal level. The com- 
puter will act as a buffer between the experimenter or the control program 
and the actual accelerator controls. For the simulation, a functional form 
for the accelerator current in terms of the signal voltage level will be 
needed, however the accelerator voltage display will be the same as the 
selected value. 

2.2.2 Electron Accelerator 

Control functions are necessary for the electron beam source, the beam 
focusing lenses and to establish the beam ejection direction. Figure 2.2 
shows a conceptual block diagram of the electron accelerator control scheme. 

The experimenter will be able to select either a voltage or a current 
for the cathode heater which provides the source of electrons for the beam. 

The voltage and current values will both be displayed and he will be able 
to select new values from the keyboard. The beam current will be determined 
by the voltage level of a control grid which can be controlled by either the 
keyboard or a preplanned control program. The accelerating potential will be 
controlled in a similar fashion. Displays for each of these source functions: 
cathode current, cathode voltage, control grid voltage, control grid current, 
and accelerating potential will be required. The experimenter can view each 
display and select a new value if necessary. Here again the computer acts 
as a buffer between the experimenter or the preplanned program and the con- 
trolled instruments. 

The focusing of the beam is done by two lenses. The first is a di- 
verging lens and the second is a converging lens. Each of these lenses is 
controlled by a voltage which is set by the experimenter. Beam positioning 
Is; accompl ished by running a current through two separate coils, one for 
movement in the Shuttle X-Z plane and the other for the Y-Z plane. The 
experimenter determines the magnitude of the current in each of the coils 
either from the keyboard, if the beam is not to be scanned rapidly, or with 
a preplanned program that is able to use the computer to impress rapid cur- 
rent variations or the coils. In these cases as nearly always the computer 
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acts as a buffer between the controlling agent and the controlled equipment. 
The computer aiso is required to record aii command settings and dispiayed 
feedbacks . 

2.2.3 The MPD Arc 

The operation of the MPD arc is controiied in two steps. First an 
experimenter determined amount of gas is iet into the arc chamber and then 
the eiectricai power bank is fired through the gas. Because the time be- 
tween these two events is very short, a prepianned program wiil be used to 
fire the MPD arc. The experimenter wiii determine gas pressures and bank 
voitage and sequence times and then modify the firing program accordingiy. 

He wili then input the program Into the computer together with a firing 
time. The computer will then send out the proper commands at the proper times 
The experimenter will need to see displays of the chamber pressure and the 
arc voltage and current. A conceptual block diagram of the MPD arc control 
scheme is given in Figure 2.3. As with the ion accelerator the beam direc- 
tion will be determined by using the Shuttle attitude system. 

2.2.4 Gimbaling System 

As stated at the beginning of this section, there will be no gimbaling 
of the accelerator system as it is presently planned. 

2.2.5 Accelerator Pointing System 

The accelerator system will be hard mounted on a pallet. The electron 
accelerator will have magnetic pointing coils as described above and the ion 
accelerator and MPD arc will be pointed by turning the whole Shuttle in the 
proper di recti on. 

2»2.6 Accelerator Deployment 

The only sort of deployment that can be considered as part of the 
accelerator system is the possible use of vacuum tight covers over the exit 
apertures of the accelerators. If such covers are used a simple experimenter 
command of open or close will be processed by the computer. Verification of 
the opening or closing would be either visual or by sensors mounted on the 
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2.2.7 Accelerator Power System 

The power systems for each of the three accelerators have been combined 
Into one energy storage bank and two power processor. The first processor 
takes 28 volt power and converts It to 500 volts. Commands and controls for 
this unit are needed to regulate the output voltage level and current. The 
experimenter will set the voltage level and current limit from the keyboard. 
The computer will keep check on the current and voltage and display their 
present values. The second power processing unit converts 500 volt power to 
power at >1,000 volts for use in the Ion and electron accelerator. Here 
again the experimenter will set the voltage and current levels and the com- 
puter will check on them and display the present values. Finally the MPD 
arc draws its power from the storage bank through a solid state switch which 
will be controlled by the computer as it runs through the preplanned MPD arc 
firing program. Figure 2.h shows a conceptual block diagram of the accelera- 
tor power system control scheme. 

2.3 HIGH POWER TRANSMITTER/ ANTENNA SYSTEM 

I While there are physically three separate power ampl if iers for the 
transmitter corresponding roughly to the three frequency ranges: 300 Hz to . 
200 kHz, 200 kHz to 2 MHz and 2 to 20 MHz, the control functions will be 
identical. 

2.3.1 Transmitter Systems 

Control of the transmitter splits into two parts, control of the trans- 
mitted wave form and control of the output power. The power levels will be 
S(pt first by choosing a maximum power consistent with the transmitter capa- 
bility and the experimental requirements, and then by means of selecting an 
attenuation value. The experimenter must choose the proper power levels; 
and the computer must activate the transmitters and attenuators, and display 
and record their outputs. Figure 2.5 shows a conceptual block diagram of 
this control scheme. 

The wave form of the transmitted waves can be generated by a separate 

! 

piece of special purpose equipment. In general the frequencies of the waves 
are too large for the spacelab computer to simulate them. However, in 
simple cases the experimenter can use the computer to set up the special 
purpose equipment. In this case he will select the carrier frequency from 
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the keyboard, select either FM or AM modulation and the frequency of the 
modulation. It may be possible to let the experimenter select other experi- 
ment specific choices from the computer keyboard. Alternatively, the ex- 
perimenter will set up the desired waveform using only the special purpose 
equipment and the computer will do no more than record all of the selected 
values. Figure 2.6 shows a conceptual block diagram of the waveform gen- 
erator control scheme In which the computer keyboard 1$ used. 

Off-on controls for all transmitter associated equipment will be done 
by the experimenter from the computer keyboard. 

2.3.2 1,000 Feet or 300 Meter Dipole Antenna 

The length of the antenna will be controlled by an extension motor. 

The experimenter will select an antenna length. The computer will signal 
the motor to start and will run the motor until the antenna length sensor 
reaches the selected value when the motor will be stopped. One further 
control function Is associated with the antenna system. This Is the matching 
of the transmitter and antenna Impedances. This function will probably be 
automatic once the experimenter has selected a carrier frequency and an 
antenna length. The computer will read the frequency and antenna length 
and select the proper matching network. Figure 2.7 shows a conceptual block 
diagram of the antenna control system. 

2 . 3.3 Transmitter Deployment System 

The transmitters and antenna are mounted to a pallet. The only de- 
ployment Is the extension of the antenna which Is treated In Section 2.3.2 
above. 

2.h BOOM A, DIAGNOSTICS BOOM 

2.4.1 Boom A, Deployment System 

The deployment of Boom A Is accomplished by having the experimenter 
select the desired boom length and turn the boom extension motor on. The 
computer will read the boom length sensor and stop the extension motor when 
the selected value is reached. Correction for length overshoot can be 
built Into the computer software. The deployment system for Boom B Is 
identical to that of Boom A and a conceptual block diagram of both systems 
is shown In Figure 2.8. This simulation will use a simple linear time de- 
pendence to display the boom length. 
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Figure 2.7 Conceptual Block Diagram of Antenna Control 
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2.4.2 Boom A Maneuvering Function 

It is necessary that Boom A be abie to maneuver in two directions. 

If when the boom is stowed Its extendable axis corresponds to the Shuttle 
Z axis, then denote the boom X and Y axes as those which correspond to the 
Shuttle X and Y axes, respectively. The two maneuvering direction corres- 
pond to angular rotations about the Shuttle Z axis and the boom axis. The 
experimenter or a preplanned program selects the desired value for each of 
these angles, and starts the rotator motors. The computer tracks the rota- 
tion sensors until the selected values are reached and stops the motors. 

In this simulation the angles are assumed to be linear functions of time 
and the computer simply displays the rotation angle as a constant multiple 
of the time that the rotating motors have been running. This simulation 
also assumes that the two angles can be varied independently of each other. 
In actual flight, It may be necessary to work out a software program which 
moves the booms in tandum because of constraints put upon their motion by 
Shuttle angular momentum considerations. 

2.4.3 Gimba led Platform 

The platform at the end of Boom A must be able to rotate about three 
separate axes. The rotations are selected by the experimenter and executed 
by the computer in the same manner as Boom A maneuvers described in Section 
2.4.2 above. The platform must also be able to remain at a fixed angle with 
respect to the Earth's magnetic field, the Shuttle velocity vector or an 
inertial frame. Therefore software must be provided so that the computer 
can maintain these fixed positions. Figure 2.9 shows a conceptual block 
diagram of the gimbaled platform control scheme. This simulation will dis- 
play the platform angles as 1 inear time functions and each angle will be 
Indepenjdentl y variable. 

2.4.4 5 Meter Sub -booms 

The only controllable function of the sub-booms are their length which 
will be controlled in the same way as the Boom A or B extension described in 
Section 2.4.1. However, there will be no need to select the length of the 
sub-booms. They will either be extended or retracted. 
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2.4.5 Electric Dipole Extension 

The control functions of the electric dipole antenna are exactly the 
same as for the 1,000 foot dipole antenna described in Section 2.3.2. 

2.4.6 Power Supply 

The power supply on Boom A requires only on-off control and displays 
of the voltage level, the charging current and the discharging current. 

2.5 BOOM B, PERTURBATION BOOM 

2.5.1 Boom B Deployment System 

The Boom B deployment system is exactly the same as the Boom A deploy- 
ment system described in Section 2,^.1. 

2.5.2 Boom B Maneuvering Function 

Boom B and Boom A have nearly identical maneuvering functions. There 
is a possibility that Boom B will be able to rotate about a third axis — 
its own Z axis — in addition to the two axes required for Boom A. Figure 
2.10 shows a conceptual block diagram of both Boom A and Boom B maneuvering 
control functions where a two axis rotation is assumed for Boom A and a three 
axis rotation for Boom B. 

2 . 5.3 Target Deployment and Retraction 

The target will be deployed by inflating it with gas. The targets 
may be made from an elastic material so that they can deflate themselves 
when a valve is opened, or they may be rolled onto a storage drum with the 
gas being ejected in the process, or finally, the targets may be jettisoned 
from the Shuttle when the experiment is finished. For this simulation we 
have assumed that target returns to its stowed position after a deflation 
gas valve is opened. The experimenter will be required to control both in- 
flation and deflation gas valves using the computer display of the target 
pressure to guide him. Figure 2.11 shows a conceptual block diagram of the 
target deployment and retraction control scheme. 

2 . 5.4 Wave Generator 

The wave generator on Boom B is used for experiments in which plasma 
waves are propogated from one boom to the other. It consists of a trans- 
mitter and an antenna, and is nearly identical to the high power trans- 
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mitter/antenna system described in Section 2.3> The output power is less 
and the antenna length is less but in terms of the control functions this 
wave generator and the equipment described in Section 2.3 are identical. 

2 . 5.5 Low Energy Electron Gun 

The low energy electron gun resembles the electron accelerator. Its 
cathode must be heated and there must be an extractor grid whose potential 
controls the electron flux leaving the gun. The experimenter will control 
the cathode heater current, and the control grid voltage. Displays of the 
cathode heater current and voltage and the grid current and voltage will 
be made by the computer. Figure 2.12 shows a conceptual block diagram of 
the low energy electron gun control scheme. There will be no beam focusing 
or magnetic bending. 

2.6 DEPLOYABLE SATELLITES 

The entire command and control system for a satellite can be extremely 
complicated expecially for a satellite of the complexity planned for the AMPS 
Shuttle payloads. Therefore, the following description applies only for the 
simple controls that were necessary for the simulation of experiments. 

2.6.1 Deployment Mechanism 

The deployment of the satellite consists, for the most part, in checking 
out the proper functioning of the satellite systems and Instruments before It 
is ejected from the Shuttle. The direction of ejection will be controlled by 
orienting the Shuttle itself and launching the satellite along the Shuttle 
plus Z-axis. Once the satellite has been checked and the Shuttle oriented, 
the experimenter will select the desired launch time and arm the ejection 
mechanism. The computer will eject the satellite at the correct time. How- 
ever, the experimenter can override the computer at any time. Figure 2.13 
shows a conceptual block diagram of the satellite deployment control scheme. 

2.6.2 Television System 

The television system requires extensive control by the experimenter. 

It will be assumed that the system Is mounted on a two axis scan platform 
so that the TV system can be targeted without changing the attitude of the 
satellite. The experimenter will select the desired angular values for the 
two axis and the command will be telemetered to the satellite scan platform 
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motors. The field of view of the TV system will be controllable and the 
magnification also. The experimenter will have on-off controls and the 
ability to change the photometric sensitivity. Finally, the wavelength of 
the received light will be selectable by changing optical filters. Figure 
2.14 shows a conceptual block diagram of the satellite television system 
control scheme. The actual television picture data will not be processed 
through the spacelab computer because the data rate is too high. 

2 . 6.3 Transponder 

The only experimenter control of the satellite transponder is to turn 
it on or off. This function can be done either through the satellite— 
Shuttle umbilical line or by means of the telemetry system. Figure 2. IS 
shows a conceptual block diagram of these two methods of turning the trans- 
ponder on or off. 

2.6.4 Telemetry System 

The only control the experimenter has on the satellite telemetry sys- 
tem will be to turn it on and off through the Shuttle-satellite umbilical 
system. Figure 2.16 shows a conceptual block diagram of the satellite 
telemetry system control scheme. 

2 . 6.5 Ranging and Control System 

The satellite range, attitude and spin rate will be controlled by 
teilemetry from the Shuttle. The experimenter selects a range or alterna- 
tively a geographic position for the satellite and the computer controls 
the satellite thrusters in order to bring the satellite to the desired po- 
sition. The attitude thrusters are then used to bring the satellite into 
the chosen orientation. Finally, the satellite is spun faster or slower to 
suit the experimenter. Figure 2.16 shows a conceptual block diagram of the 
satellite ranging and control system. 

2.7 DEPLOYABLE UNITS 

The control of the deployable units is divided Into two parts. The 
first is concerned with the ejection of the unit and is exactly like the 
procedures described in Section 2.6.1. The second facet of the control of 
the deployable units involves the telemetering of commands from the experi- 
menter to the unit. A conceptual black diagram of the deployable unit con- 
trol shceme is shown in Figure 2.17. 
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Figure 2.15 Conceptual Block Diagram of Satellite Transponder Control 
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Figure 2.16 Conceptual Block Diagram of Satellite Ranging and Control System Control 












Figure 2.17 Block Diagram of Deployable Unit Control 
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3.0 EXPERIMENT SIMULATION 


A complete description of the simulation requirements of six experi- 
ments are presented in the following sections. 


Experiment 1 
Experiment 2 
Experiment 3 
Experiment k 
Experiment 5 
Experiment 6 


Electromagnetic Wave Transmission 

Passive Observation of Ambient Plasmas 

Ionospheric Measurements with a Subsatellite 

Electron Accelerator Beam Measurements 

Lidar Trace of Acoustic Gravity Waves in the Sodium Layer 

Determination of the Wake of a Test Body 


For each experiment simulation, the responses of a subset of sensors and 
controllable systems listed in Table 1.0 and 2.0, respectively are described. 
The sensor and detector response are presented as a logical function of the 
relationship between the detector orientation with respect to its velocity 
vector, the orientation and velocity vector of the Shuttle, and its location 
in orbit. This includes the following factors: 

a. The pointing, control, data acquisition and display of the remote 
sensing platform and the detectors and sensors thereon. 

b. The deployment, pointing, operation and display of the gimbaled 
accelerator systems. 

c. The deployment orientation and operation of the high power trans- 
mitter/antenna system. 

d. The deployment, maneuvering, operation, control and display of 
Booms A and B including the operation and display of all instru- 
ments, sensors, detectors, and associated subsystems on each. 

e. The deployment, tracking, and display of all deployable units. 

f. The deployment, tracking, control , operations , and display of 
the two subsatellite units. 

This section is subdivided into six parts, one for each of the experi- 
ments. Each part in turn is divided into three subparts: An experiment 

description, an experiment procedure and a theory and background section. 

The data formating requirements for each experiment are given in the 
appendix. 
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3.1 SIMULATION OF THE ELECTROMAGNETIC WAVE TRANSMISSION EXPERIMENT 

3 , 1.1 Experiment Description 

Experiment description uses the outputs of the following Instruments 
described in Section 1.0. 

1. One meter loop antenna 

2. Short electric dipole antenna 

3 . Triaxial fluxgate magnetometer 

4. 33 meter electric dipole antenna 

5 . Cylindrical electron probe 

Furthermore in response to the requirements of Task 2 the experiment 
exercises the following controllable systems: 

1. The high power transmitter antenna system 

2. the Boom B perturbation boom 

3 . The Boom A diagnostics boom with platform. 

The outputs of the instruments and controlled systems are flow charted In 
Section 4.0 as required by Task 3 and Integrated Into the experiment simula- 
tion as required by Task 4. 

This experiment simulation has been carried out in great detail in 
order to demonstrate the flexibility and control that can be provided m a 
CVT simulation. This detail is particularly evident in the theory and in 
the euqation derivations. Furthermore each display presentation of instru- 
ment operation has been backed up with another display showing computer de- 
fault values for experimenter choices, thus demonstrating the systems capa- 
bility to handle default values. 

A wave transmission experiment would be performed early In an AMPS 
mission. The theory of 1 inear electromagnetic wave transmission is well 
understood so that an experiment like this one would give experimenters 
early practice In the use of the active experimental capability of the AMPS 
payload. This particular simulation experiment provides a test of the wave 
transmission and receiving equipment and could be carried out in a variety 
of orbital situations to ascertain that the equipment is functioning properly. 

The simulation experiment consists of transmitting a wave of known fre- 
quency and amplitude from one center-fed dipole to a short dipole antenna and 
a loop antenna which measure the electric and magnetic fields of the trans- 
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mitted wave. The transmitter is stepped in frequency from a frequency 
several times larger than the local plasma frequency down to a value very 
close to the plasma frequency. The experimental result is the ascertainment 
of the wave vector, k, as a function of transmitter frequency and ionospheric 
plasma frequency. 

During the wave transmission experiment the cylindrical electron probe 
and the fluxgate magnetometer are making measurements of the ionospheric 
electron density and temperature, and the Earth‘s magnetic field. Because 
the transmitting and receiving antennas are mounted on separate booms, the 
experimenter has the ability to vary the distance and angle between them. 

A sample experimental procedure is presented in order to introduce the ex- 
perimenter to the simulation format. Once he has learned how to get results, 
he will be able to vary some parameters at will. The simulation is limited 
in that the following quantities cannot be varied: 

1. The electron density and temperature 

2. The magnetic field value 

3. The Shuttle orbital parameters. 



However, there still remains a wide variety of experimental conditions avail- 
able to the experimenters. 

3*1*2 Experiment Procedure: Wave Transmission Experiment 

EPl.l Hold Shuttle in orientation with Y-axis along ambient magnetic 
field direction (orbiter controT system) 

EP1.2 Set boom B 

EPl.2.1 Set length to 50 meters relative to Shuttle 
EPl.2.2 Set 0-angle at 30° relative to Shuttle 
EPl.2.3 Set 0-angle at 0° relative to Shuttle 

EPl . 3 Set boom A 

EPl.3.1 Set length to 50 meters relative to Shuttle 
EPl. 3. 2 Set 0-angle at 30° relative to Shuttle 
EPl. 3. 3 Set 0-angle at 180° relative to Shuttle 
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EPI.4 Set boom A gimbaled platform 

EPI.4,1 Platform in magnetic autolock mode 

EP1.4,l.l Set short electric dipole antenna axis 

parallel to ambient magnetic field [must 
display orientation] Numerical Table 

EPI.4.1.2 Set autolock [Computer must check magnetic 
orientation and control] 

EPI.5 Set wavegenerator length at 6 meters tip-to-tip 

EPI.6 Generate background plasma data 

EPI.6,1 Turn Boom A Power Supply ON 
EPI.6.2 Cylindrical electron probe 

EPI.6.2.1 Turn instrument ON 

EPI.6.2.2 Check sweep voltage waveform CRT display 

EPI.6.2.3 Check current versus voltage CRT display 

EPI.6.2,4 Check probe orientation display probe not 
parallel velocity vector 

EPI.6,2.5 Perform electron temperature calculation 
EPl.6.2.5.1 Input model 
EPl ,6.2. 5. 1 . I Display Te = l/lOeV 
EPI.6.2.5.1.2 Record 

EPl. 6. 2. 6 Perform electron density calculation 

EPl. 6. 2. 6. I Display n^ ■ 3x10^ electron/cm3 
EPl. 6. 2. 6. 2 Record 

EPl. 6. 3 Fluxgate magnetometer 

EPl. 6. 3.1 Turn fluxgate ON 

EPl. 6. 3. 2 Display magnetic field values 

EPl. 6. 3. 3 Calculate electron cyclotron frequency 

EPl. 6.3. 3.1 Display cyclotron frequency #F ■ 
840 kHz 

EPl .6. 3. 3. 2 Record 

EPl. 6. 3. 4 Record magnetic field values 

EPl. 6. 4 Calculate plasma frequency 

EPl. 6. 4,1 input electron density 3x10^ 

EPl. 6. 4. 2 Display plasma frequency F^ » 4.9 MHz 
EPl, 6.4. 3 Record 
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EP1.7 Generate electric field background survey from short electric dipole 

EPl.7.1 Display short electric dipole is parallel to magnetic field vector 
EPl.7.2 Connect electrjc field receiver to short dipole antenna 
EPl.7.3 Set receiver bandpass 

EPl.7.3.1 Upper level 20.5 MHz 
EPl.7.3.2 Lower level 3.9 MHz 
EPl. 7 . 3.3 Square rolloffs 
EPl.7.^ CRT display receiver output 

EPl. 7.5 Adjust gain so that background level exists >0 at all 
frequencies 

EPl. 7.6 Record electric background survey 

EPI.8 Generate magnetic field background survey from loop antenna 

EPl.8.1 Check loop antenna perpendicular to magnetic field vector 

EPl.8.2 Connect magnetic field receiver to loop antenna 

EPl. 8.3 Set receiver bandpass 

EPI.8.3.1 Upper level 20.5 MHz 

EPl. 8 . 3.2 Lower level 3.9 MHz 

EPl. 8 . 3.3 Square rolloffs 

EP1.8.A CRT display magnetic receiver output 

EPl.8.5 Adjust gain so that background level exists >0 at all 
frequencies 

EPl.8.6 Record magnetic background survey 

EPI .9 If background levels are acceptable continue, if not, stop 

EPl.lO Activate transmitter 

EPl.10.1 Connect 33 meter antenna to transmitter 
EPl.lO. 2 Turn transmitter to Ready 
EPl.10.3 Set attenuation to maximum 
EPl.lO. A Set carrier waveform to sine wave 
EPl.lO. 5 Set carrier frequency to 19.60 MHz 
EPl.10.6 Turn transmitter ON 

EPl.ll Receivers 

EPl.11.1 Set electric field receiver bandpass at 19.61 and 19.59 
EPl.ll. 2 Set magnetic field receiver bandpass at 19.61 and 19.59 
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EP1.12 Adjust transmitter attenuation so that received signal level Is 
100 times background level 

EPl.12.1 Display received Electric Field Signal 

EPl.12.2 Adjust attenuator 

EP 1.12.3 Record all levels and settings 

EP1.13 Measure fields 

EPl.l3«l Electric field receiver 

EP1.13.1.1 Digitize mean square electric field 
EP1.13>1.2 Generate data for 30 seconds 
EPI.I 3 .I .3 Transmitter to Ready 
EPl. 13 . 1 .^ Generate mean square electric field 
EPl. 13 . 1.5 Display mean square E-field 

EP1.13*2 Magnetic field receiver 

EPl. 13«2.1 Digitize mean square magnetic field 
EP1.13>2.2 Transmitter ON for 30 seconds 
EPl. 13 * 2.3 Transmitter to Ready 

8 EP1.13.2,it Generate mean square magnetic field 

EPl. 13^2. 5 Display mean square B-field 

EPl. 14 Calculate k- value 

EPl. 14.1 Call for special program "calculate k-value" 

EPl. 14. 2 Display k value 

EPl. 14. 3 Record all calculations 

EPl . 15 Repeat steps EPl. 10. 3 to EPl.14.7, inclusive, and step EPl. 6. 4 
frequencies as below 

EPl. 15.1 Set receiver bandwidth to F + .01 and F - .01 of transmitter 
frequency F 

EPl. 15.2 Set transmitter frequency F at 
EPl. 15. 2.1 15.5 MHz 

EPl. 15. 2. 2 10.9MHz 

EPl. 15.2.3 8.95MHz 

EPl. 15.2.4 8.0 MHz 

EPl. 15 . 2. 5 7.5 MHz 

EPl. 15 . 2.6 7.0 MHz . 

I EPl . 15 . 2.7 6.75 MHz 

EPl. 15.2. 8 6.50 MHz 
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EPl.15.2.9 6.25 MHz 

EPI.15.2.10 6.00 MHz 
EPl.15.2.11 5.9 MHz 

EPI.I5.2.12 5.8 MHz i 

EPl.15.2.13 5.7 MHz 
EPl.15.2.14 5.6 MHz 

EP1.15.2.15 5.5 MHz 
EPl. 15.2.16 5.^ MHz 

EPl.15.2.17 5.3 MHz 
EPl. 15. 2. 18 5.25 MHz 

EPl. 15. 2. 19 5.20 MHz 

EPl. 15.2.20 5.15 MHz 

EPl. 15. 2. 21 5.10 MHz 

EPl. 15.2.22 5.05 MHz 

EPl. 15.2.23 5.00 MHz 

EPl. 15. 2. 24 4.95 MHz 
EPl. 15. 2. 25 4.90 MHz 
EPl. 15. 2. 26 4.85 MHz 
EPl. 15. 2. 27 4.8 MHz 
EPl. 15. 2. 28 4.7 MHz 
EPl. 15. 2. 29 4.6 MHz 

EPl. 15. 2. 30 4.5 MHz 
EPl. 15. 2.31 4.0 MHz 

EPl. 16 Terminate Experiment 

EPl. 16.1 Transmitter OFF 
EPl. 16. 2 Electric Receiver OFF 
EPl. 16,3 Magnetic Receiver OFF 
EPl. 16.4 Disconnect Antennas 

EPl. 16. 4.1 Transmitter from 33“meter Antenna 
EPl. 16. 4. 2 Receiver from Electric Dipole 
EPl. 16. 4. 3 Receiver from Loop Antenna 
EPl. 16.5 Reel in 33-meter antenna 
EPl.16.6 Turn Fluxgate Magnetometer OFF 
EPl. 16.7 Turn Cyl indrical Probe OFF 
EPl. 16. 8 Stow both larger booms. 
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3-1-3 Theory and Background for Simulation of the Electromagnetic 
Wave Transmission Experiment 

3-1-3-1 AMPS instruments and Equipment Required . The following AMPS 
instruments and equipment are used in this experiment: 

1) Transmitter 

2) Wave generator (variable length antenna) on Boom B 

3) Electric field receiver with dipole antenna on Boom A 

4) Magnetic field receiver with loop antenna on Boom A 

5) Boom A with glmbaled platform 

6) Boom B 

7) Fluxgate magnetometer on Boom A ■ 

8) Cylindrical probe on Boom A. 

3. 1.3. 2 Simulation Controls. In the experiment simulation and experi- 
meter can control the following parameters over the given ranges: 


1) 

ONrOFF functions for every instrument 

ON-READY-OFF 

2) 

Transmitter frequency 

3 to 25 MHz 

3) 

Transmitter amplitude and/or power 

0.1 to 1000 watts 

4) 

Background E-M noise level 

10 to 10’^ volts 

5) 

Wave generator (transmitter) antenna 
length 

0-33 meters 

6) 

Electric field receiver bandpass 

Any at frequencies 
from 3 to 25 MHz 

7) 

Electric field receiver gain 

0 to 80 dB 

8) 

Magnetic field receiver bandpass 

Any frequency 
from 3 to 25 MHz 

9) 

Magnetic field receiver gain 

0 to 80 dB 

10) 

Boom A length 

2-50 meters 

11) 

Boom A ® angle 

0-90° 
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0-360° 

0-360° 

0-90° 

0-360° 

2-50 meters 

0-90° 

0-360° 

3. 1.3*3 Simulation Displays . The following simulation displays are 
provided and are shown in Section 3*1 -3.^. 

1) Wave transmission experiment display 

2) Boom B set-up display 

3) Boom B set-up display default 

4 ». 4) Boom A set-up display 

5) Boom A set-up display default 

6) Platform set-up display 

7) Platform set-up display default 

8) Cylindrical electron probe set-up display 

9) Cylindrical electron probe set-up display default 

10) Cylindrical electron probe voltage vs time 

111) Cylindrical electron probe current vs voltage 

112) Fluxgate magnetometer set-up display 

13) Fluxgate magnetometer set-up display default 

li4) Electric field receiver set-up display 

lib) Electric field receiver set-up display default 

16) Magnetic field receiver set-up display 

17) Magnetic field receiver set-up display default 




12) Boom A (}> angle 

13) Platform x angle 

14) Platform ^ angle 

15) Platform n angle 

16) Boom B length 

17) Boom B ® angle 

18) Boom B 4> angle 
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18) Transmitter set-up display 

19) Transmitter set-up display default 

20) Wave transmission experiment set-up display 

211) Wave transmission experiment E set-up display default 

22) Wave transmission experiment set-up display 

23) Wave transmission experiment B set-up display default 

24) E data acquisition display 

25) E data acquisition display default 

26) B data acquisition display 

2 

27) B data acquisition display default 




3.1.3-4 Special Purpose Equipment Required for Simulation . In order 
to perform this simulation the following special purpose equipment must be 
added : 


1) Moderately fast, oscilloscope, 200 MHz ability, preferably 
a storage type 

2) A sine wave source able to cover the frequency range 
from 3 to 25 MHz at 10 volts output into 1 Mq. 

3) A noise source able to output in the frequency range 
3 to 25 MHz with voltage amplitudes on the order of 
0.01 to 0.1 volt into 1 MJ2. This source need not be 
genuinely random, but it should at least look noisy. 

4) A wave frequency spectrum analyzer having a bandwidth 
resolution on the order of 10 to 100 kHz in the fre- 
quency range from 3 to 25 MHz. The analyzer must have 
its own display or be compatible with the oscilloscope. 

5) A switching panel able to be controlled by the computer 
and used to connect the noise or sine source or spectrum 
analyzer outputs to the oscilloscope or to the spectrum 
analyzer. 

6) the noise source and the sine wave generator must be 
capable of having their amplitudes controlled by the 
computer. The sine wave generator frequency must also 
be computer controlled. Some piece of hardware will 
probably be needed to interface these sources with the 
computer. 
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3. 1.3.5 Assumptions Used in Simulation . In order to do the simulation 
of the wave transmission experiment certain assumptions must be made con- 
cerning the Shuttle orbit, the state of the ionosphere and other determining 
factors. The following is a list of these assumptions for this one 
experiment: 

1) The earth is a perfect sphere of radius 6371 km. 

2) The earth's magnetic field is the result of an earth- 
centered dipole pointed at the north geographic pole 
and having a dipole moment of 

8.07 X 10^^ gauss-cm'^ = 8.07 x 10 ^ teslas-m . 

3) The Shuttle is in a circular orbit of zero inclination 
at an altitude of 400 km above the earth. 


4) The ionosphere is uniform at this altitude throughout 
the entire orbit. 

5) The ionospheric electron density at this altitude is 
3 X loil electrons per m^. 

6) The ionospheric electron temperature at this altitude 
is 1161°K or 0.1 eV. 

7) The level of background electromagnetic radiation is 
whatever the experimenter sets by using the noise source. 

8) There are no other background, interference or environ- 
mental functions that affect this experiment, or that 
need be specified. 

9) The phases of all oscillations are neglected. Assump- 
tions made in boom simulation were: 

a) A motor running an angular variable runs at the 
rate of one degree per second in realtime. 

b) A motor running a length variable runs at the rate 
of 0.5 meter per second in realtime. 

c) Lengths are known to 0.1 meter and angles are known 
to 1.0 degree. 

3. 1.3.6 Wave Transmission Experiment Theory . Simplifying assumptions 

1) Everywhere in space the electric and magnetic fields are 
as they would be in a vacuum except that the plasma 
effects are accounted for by setting 
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w _ 


^ /I - (Wp/u))^ 


Instead of the usual 


u)/k = c and Wp = plasma frequency 


/4,m 


ectron density 


2) The length of the transmitting antenna is small compared 
to either the wavelength or the separation between the 
transmitting antenna and the receiving antenna. 

3) The transmitting antenna is treated as a simple center 
fed dipole. 

4) Under the conditions of the experiment u ^ , w are small 
and the dispersion relation for both the ordinary and 
extraordinary waye modes is given by 


o)/k = c/(l - {wp/co)^j 


2 \ 1/2 


when CD > ojp, and no wave exists for cd < cDp. 

In gaussian units the dipole moment of the transmitting antenna is 
given by 

i I. 


P = 


od - 
2 (d ®y 


where 


d = antenna length 

1 = AT 

CD = driving frequency 

I = current leaving the transmitter 

0 

e = unit vector along the antenna 
^ length 


( 1 ) 


or 


p = p e 
K ►'o y 
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The figure shows the geometry of the experiment: 



RECEIVING 

ANTENNA 


DJ SHUTTLE 


The coordinate system used in the theory is centered on the transmit- 
ting antenna (the tip of Boom B) and the antenna is pointed along the 
y-axis of this coordinate system. 

The electric (E) and magnetic (B) fields at a point x, y, z are 
given by 

. 2 , . 2 , 


E = P^e 




( 2 ) 


and 


i(kr-o)t) 
c r V 


- ik)(-zg^ + 


( 3 ) 


where i = k = wave number; c = speed of light in vacuum. 


r = /x^+y^+z^ = distance from center of transmitting antenna 
to point of measurement 


i I 


p^ = dipole moment = -s— 

'^0 Z(JL) 


od 


u = wave frequency 
= transmitter frequency 


®Y» ®\/» are unit vectors. 

^ V ^ 
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The electric field value actually measured by a short dipole at point 
(x,y,z) is given by 

E j = Real(E • R) where R is a unit vector pointed 

along the receiving dipole antenna 


The magnetic field value actually measured by a loop antenna is given by 


/V /V 

B . = Real(B • L) where L is a unit vector normal to 

the antenna plane 


On the Shuttle the antenna orientation is such that L • R = 0. 

Thus we now know what E and B fields will be measured at any point in 
space when the transmitting antenna is being driven with a current Iq at a 
frequency w. 

In doing the experiment to be simulated one wishes to determine what 
the wave number of an ordinary electromagnetic wave is as the frequency of 
the transmitted wave approaches the plasma frequency from the high side. 


This determination is most easily done if one looks at the ratio of 
the time averaged square of the electric and magnetic fields measured by 
the receivers: 


"^^measured^ 

<B^ .> 

measured 


2 2., 4 4, 


> _ c^(9+3k^r^+k^r^). 


w^r^(Uk^r^) 


2/2 ..2, 


C0S^Y£ 



where we have assumed that the loop and dipole receiving antennas are 

both located at point (x,y,z) with respect to the center of the transmitting 
2 2 2 2 

antenna and r = x + y + z ; c = speed of light in vacuum; to = transmitter 
frequency; and k = wave number. 

The wave number, k, is given by 

k = 

where 



6 8 

o)p = plasma frequency = 30.77 x 10 rad/sec and c = 3 x 10 meters/sec. 
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Now at any given transmitter frequency (o), the relative position 
between the transmitting and receiving antennas (x,y,z), the transmitting 
antenna tip-to-tip length (d), the transmitter current (Iq)» the receiving 
antenna orientation (R), and the electric and magnetic signal strengths 
observed by the receivers can be calculated and displayed by controlling 
the amplitude of the sine wave generator. 


The transmitter radiated power is given by 
3 k^I ^d^ 

’’cgs ' ' '"rad * "ant>^< 


Total Power 




+ 1 .955 X lO"'-^ 

I CLii 


Rant^6m tip-to-tip) = 0.1 76n 
in cgs units 


( 8 ) 


where it is assumed that the antenna is only resistively coupled to the 
environment. 

3 . 1 . 3 . 7 Angular Transformations for Boom Movement Used in Simulation . 
The following assumptions about the movement of the booms are made: 

1) The booms are able to extend in length from 2 to 50 meters. 

2) They are able to rotate about the Shuttle's Z axis by an 
angle ({>. 

3) They can rotate about their own y-axis, which is defined 
as the Shuttle Y-axis when = 0, by an angle®. 

Suppose the Boom B tip has Shuttle coordinates Lg, Og, (|>g. Further- 
more, suppose there is another coordinate system frozen into the tip of the 
boom such that when Lg = = 0 that axes x,y,z are aligned with the 

Shuttle X,Y,Z axes then the boom tip has Shuttle coordinates 

4 “ ’"B ^’’'Pb *’"4 
4 ' 4 '^pppb . 
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Boom A has coordinates 


sine^ cos<^^ 

sine^ stn.^^ (10) 

h ' '■A “®®A 

In the theory we have used the quantities x,y,z to denote the coordinates 
with respect to the Boom B tip of the receiving antennas at the tip of 
Boom A. Coordinates are determined in the following manner. First, trans- 
form from the Shuttle coordinates to the x,y,z Boom B coordinates. The 
rotational transformation matrix to carry the Shuttle into Boom B is 


/cOS6g COS4>g 

cos 0 g sin(j)g 

-Sln6g\ 


-sin(j)g 

I 

COS((.g 


(11) 

\sin6g cos<j)g 

sin<j)g sin0g 

coseg / 



and the quantities used in the theory are then found from 



/ 

y 

= 

\zl 



cos0g sitting 
COS<j)g 


-sin6g\ /l^ sine^ coS(^\ 


0 

cose, 


sine^ sin 4 >^ 


ysinOg cosiJ>g sin<j)g sineg cosQ^J 




l°\ 


A ""®A 


f 

/ 


(12) 


X = +Lg sine^ cos4>^ cos6g cos(|)g + sine^^ sin^^ cos9g sin4»g - cose^ sin8g 
y = -L^ sine^ sincj)^ sin«{)g + sine^ sinc^^ cos(j>g 


(13) 


z = -Lg + sine^ cos(|)y^ sinOg cos(j>g + sine^ sin^^ sTn(j>g sinOg 


+ cose^ cos 0 g 
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The distance from the transmitter to receiving antenna is found from 

r = /x^+y^+z^ 

One further problem is that measured signals depend upon the various angular 
orientation by Equations 4 and 5. 

The unit vectors R and L which are expressed in a coordinate system 
fixed with respect to the gimbaled platform on Boom A must be transformed 
to boom coordinates. 


THI xj, X^, xj OF »OOM A AND X, Y, Z OF »OOM B 
CO«MSFONO IWNTICALLY TO THE SHUnU X, Y, Z 
AXIS WHIN THI BOOMS A*l IN T< STOWID POSITION. 



X’,, X'j, X-j ARE FIXED ON TIP OF BOOM A AND 
X,, Xj, Xj are FIXED ON PLATFORM 
WHEN O - X - - 0 THEY ARE IDENTICAL 

CYLINDRICAL PROBE AXIS IS Xj AND THE aUXGATE 
HAS AXES CORRESPONDING TO X,, Xj, Xj. 

R - Xj and L • Xj 

Let X,q,i|; be the usual Eulerian angles of the platform with respect to a 
coordinate system x^.x^.x^ fixed in the tip of the Boom A similar to the one 
on the B Boom. We want to find orientation of x^.x^.x^ terms of the 
x,y,z-Boom B coordinates 
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^3/. 


in Boom 
B coordinates 


Tp^ = transformation from platform to Boom A 


"2 


(14) 


itostji cosX - cosfl sinX sin\|/ -sinij; cosX - cosn sinX cos4» sinP sinx\ 


cosij^ sinX + cosf^ cosX sinij; -simjj sinX + cosn cosX cosij) -sinn cos X 
^infl sinij^ sinfi cos(i» cosf 2 


T. = transformation from Boom A coordinates into Shuttle coordinates 


(15) 


/ 




sin0« cos 


jcosQp^ cos(^f^ ‘Sin(j>^ 

cose^ sint))^ cos^p^ sin<j>^ sine^ 

-sine« 0 cose 


^sine^ 


(16) 


/ 


and Te n is already given in Equation 11. In particular, we want to find 


E»R and B*L where R = Xg» L = x^. 
^measured " Rea^E-R) = E 




(17) 


* ^z^s->b''’a^s’’^p+aIL 


where [x, !y, iz means the scalar magnitude of the x, y and z components of 
the vector. 

'0\j 

= cos0g cos<|>g|sin<j)^(sinij) SinX - cosfl cosX cosij>) 




+ sine^ cos(j)^ sinJJ cosiji 


(18) 


cose^ cos(j)^(sini|; cosX + cosq sinX 


+ cos 0 g sin(|)gjcos({>^(cosn cosX cosi|) - siruj; sInX) 
+ sincj)^ sine^ sinfi cosip 

- cose^ sin(j>^(sirnj< cosX + cosn sInX cosi|<)j 

- sin0gjsin0y^(siru|^ cosX + cosfl sinX cosip) 


+ cos0^ Sinn cosip 


] 


T (f 


= cosfgjsinct)^ sin0^ sinn cosijj 
+ cos<j)^(cosn cosX cos^; - sini^i sinX) 

- cos0y^ sin<i)^(sinijj cosX + cosn sinX cosij;)j 

- sin4ig|-cos0^ cos(})y^(siniiJ cosX + cosn sinX costjj) 
+ sin(()^{-cosn cosX costj; + sinij< sinX) 

+. sin0^ cos<|)^ Sinn cosijij 


( 18 ) 



= sinOg cos(|)gj^-cos0^ cos(}!y^{sintiJ cosX + cosn sinX cost|i) 
+ sin(j)^(sinip sinX - cosn cosX cosij;) 


+ sin0^ cos4>^ 

+ sin<{)g sin 0 g|j-cos 0 ^ sin<})^(sinij.' cosX + cosn sinX cos»|j) 
+ cos<|)y^(cosn cosX cosij/ - sin\i» sinX) 

+ sin(j)^ 

+ COS0g 
+ COS0^ 


sine^(sinii) cosX + cosn sinX cosii*) 
sinn cosij)! 


Sinn COS4; 


] 


I 


I 


I 
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3. 1.3. 8 Equations for Wave Transmission Experiment Simulation Displays 

3. 1.3. 8.1 Received Signals 

2 

3. 1.3. 8. 1.1 General Display of Received Signal Components <E > and 
<B^> . The values of mean square electric and magnetic fields displayed in 
graph form are calculated for values of all variables taken at t - 0.5 sec- 
onds. The point is plotted at t = t seconds. The value is calculated from 


<E^> " <eE 


measured’ ' sin(kr-mt+«)> = ^ 


where 6 is the phase angle and 


c2 _ 81 X 1 Power in watts ^ • ■ v . \ 

4(9 X 10^ .Uejui^r^ 


f£= Vb0 - Vb(j 


(2 - 


' 0 \| 12 



where 



are given in Equation 18 and all dimensions are measured in MKS units. 

I ‘k 


6b 



■i'J f 


“ ‘®Lsured" " <B^(s1n(kr-a.t« ')) > = ^ 

,2 _ 2.5 X 10’^® (Power in watts) f rn n t. a i i y a 


( 


3 ^ 

9 k'^d 

9x 10^ M/ 


.ue) 


»n) 




I^T 




fVb‘° 



X 




(20) 


and 




z or X 


is calculated in some fashion as was 

fO 




z or X 


in Equation 18 from Equations 11, 14, 15, and 16. 

8. 1 .2 Received Signal at Fixed EPl.O Boom Positions 

In EPl.O the received signals are measured with the boom fixed such 
that = Lg = 50 meters, 

6^ ” 6g “ 30 , <|)^ “ 180 , <j>g “0 

and 

X = ^^ = i() = 0 

In tl ifs case the mean square electric and magnetic fields reduce to 


<e2> = <e2 > 

measured 


JO 


= 2.33 X 10 *^ (P) (1 - 2 . 5 x 10 ^ + 6 . 25 x 10 °) 


^ 2.7 X 10 ^° ^ + . 176 ^ 

where P is the power in watts as chosen by keyboard 


k = 


_ /u )2 - 9.4668 X 10^4 


3 X 10 


W 


per meter 


/ vol tsV 
\metery 


0) 


= 2TTf 


where f is the transmitted frequency and 


2 2 

<B > = <B^ .> 

measured 


= 5.4 X 10 '^’ (1 + 2 . 5 x 10 ^) 




7x10’“^. 

0) 


. 176 ^ 


teslas^ 


3.1. 3-8. 2 Cylindrical Probe Angles . For experiment EPl .0 we assume 

-¥■ '' 

that the Shuttle Is moving in the minus X direction, e.g., v = X. 


The angle Q between the cylindrical probe axis (= Xg of platform) and 
V is given by 


cos Q = X . Xg = -VsVAiJ, 


= cos cos (|)^ (+sini|) cosX + cosn sinX cosiji) 
sin4>y^ (sim|7 sinX - cosn cosX cosij») 


^ sine^ cos<|)^ sinn cosip 
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3.1. 3.8.3 Fluxqate Magnetometer - Magnetic Field Components . To 
find the measured components of the magnetic field B^^ called for 

in the fluxgate magnetometer display the total B field vector, which dur- 
ing this experiment lies in the plus Y Shuttle axis, e.g., B = Bj Y, 
is transformed into the platform coordinates, x^.XgtXg* so that the field 
values at the magnetometer are given by 


h " ^|»A 




1 


®Y " 




1 -1 I 

1 -1 1 

-aVs(’KL 


( 22 ) 


(23) 


(24) 


■ n 

- 

VTa*’' 


BTOT = = 25997 GAMMA 


during the whole experiment, and 




'cose^ cos<j)^ 


-sin<j). 


.-1 ^ 
p>A ” 


^ine^ cos<)>^ 

ycosi|) cosX - cosj^ sinX sinij; 
-sinij( cosX - cosJi sinX cosij) 


V 


nn sinX 


\ 


cose^ sincj)^ -sme^ 


COS(j), 


sin(j)^ sin9^ cose^y 


cos(|j sinX + coss^ cosX si mj) sins sinij^ 


,\ 


-sim|) sinX + cosfi cosX cosij) sinQ cosij^ 
-sinJi cosX coss^ 


/ 


I 
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3. 1.3. 8.4 Angles Between Magnetic field and the Antennas or Booms . 
The angle between receiving short dipole antenna and magnetid field RB is 
I given by 

cos RB = Y • = W|hA0 

^ 'Y component 

! 

between the magnetic field and the loop antenna LB 








COS LB = Y • X3 = 


0 


Y component 


between the transmitting antenna and the magnetic field TB 


cos TB = T 


S.0 


Y component 


The cosine of the angles between the magnetic field and the coordinate 
cirections of the boom A, boom B and platform coordinate systems are given 
below: 


BOOM A 

PLATFORM 

BOOM B 

D i recti on 

Cos of Angle 

Di recti on 

Cos of Angle 

Direction 

Cos of Angle 

Ai 

Vs 1 

u Y 

A. 

Vs'^p A ° 

u Y 

X 

Tsleg 

u Y 

r— CM 
< X 

0 

VsJ 

u Y 

A 

X 2 

0 

Vs^p A 1 

u Y 

Y 

-1 0 

VbJ 
u Y 

4 

0 

Vs? 

' Y 

■ ' A 

>- 

0 0 ■— 

S' 

^ 

Z 

-1 ° 
Vb? 

Y 



3. 1.3. 8. 5 Angles Between the Various Antennas . The angle between 
the transmitting antenna and the receiving short dipole antenna TR 


cos TR = 

^ ■ Y component 

between the transmitting antenna and the loop antenna TL 

cos TL = VbVs^P^a0 

W component 

3. 1.3. 8. 6 The Actual Transmitter Power Output 

- Attenuation dB 

’’actual = (’’maxXSignal Level 10 
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3.2 EP2.0 PASSIVE OBSERVATION OF THE AMBIENT PLASMA 

3.2.1 Description 

In order to carry out active experiments in the ionosphere it is 
necessary to have an accurate knowledge of the ambient state of the ion- 
osphere. Passive ionospheric measurements have been carried out for 
nearly 40 years and during the last 20 years there have been a large 
number of in-situ satellite and rocket measurements of the ionosphere. 

The AMPS Shuttle payload must be equipped with a large set of instruments 
capable of making the necessary ambient ionospheric measurements. 

In this simulation a simple model of the ionospheric region from 
120 to 400 km is developed. The model shows local time, latitude, and 
altitude dependencies of the ionospheric state function. The outputs of 
the various ionospheric instruments on the Shuttle respond to the varia- 
tions in the model. The experimenter is able to vary the Shuttle orbit 
altitude, inclination, attitude, and the boom orientation. For this 
simulation it has been assumed that the presence of the Shuttle in no 
way affects the ability of the instruments to measure the ambient. During 
an actual flight there will exist a considerable wake region around the 
Shuttle where it will be impossible to make ambient measurements. 





■ ' : J , I " 1 I ■ ■ J ^ - I 
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3.2.2 Experiment Procedure - EP 2.0 Passive Observation of the 
Ambient Plasma 

EP2.1 Determine Shuttle Orbit 

2.1.1 Select manual control 

2.1.2 Shuttle altitude equal 300 km 

2.1.3 Set Shuttle inclination equal 57 'deg 

2.1.4 Set local time of southward moving node equal 
0000:00 

2.1.5 Set universal time of descending node equal 0000:00 

EP2.2 Determine Shuttle Attitude 

2.2.1 Select manual control 

2.2.2 Set r = 0-deg 

2.2.3 Set A = 0 deg 

2.2.4 Set A = 0 deg 

2.2.5 Hold Shuttle at these angles 

EP 2.3 Set Boom A 

2.3.1 Set length to 50 meters 

2.3.2 Set 9=0 deg 

2.3.3 Set c|) = 0 deg 

EP2.4 Set Platform 

2.4.1 Set X = 0 cleg 

2.4.2 Set ij; = 0 deg 

2.4.3 Set = 0 deg 

EP2.7 Fluxgate Magnetometer 

2.7.1 Fluxgate to GO status 

2. 7. 1.1 Fluxgate to ON 

2. 7. 1.2 Wait for Fluxgate READY 

2. 7. 1.3 Fluxgate to 60 

2.7.2 Observe magnetic field values 

2.7.3 Record all future data 

EP2.8 Rubidium Magnetometer 

2.8.1 Rubidium to GO status 

2. 8. 1.1 Rubidium to ON 

2. 8. 1.2 Wait for Rubidium READY 

2. 8. 1.3 Rubidium to GO 

2.8.2 Observe magnetic field value 

2.8.3 Record all future data 






I 


I 


I 




^ V 


: rt I'.' 
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EP2.9 Cylindrical Electron Probe 

2.9.1 Cylindrical probe to GO status 

2.9.1 .1 Probe to ON 

2. 9. 1.2 Wait for probe READY 

2.9.1 .3 Probe to 60 

2.9.2 Observe probe output 

2. 9. 2.1 Electron temperature 

2.9. 2. 2 Electron density 

2.9. 2. 3 Probe voltage versus time 

2. 9. 2.4 Probe current versus voltage 

2.9.3 Check that angle between Shuttle velocitv y^ftor 
and probe axis is greater than 10 deg 

2.9.4 Record all future data 

EP 2.10 Spherical Ion Probe 

2.10.1 Spherical probe to GO status 

2.10.1.1 Probe ON 

2.10.1.2 Wait for probe READY 

2.10.1.3 Probe to 60 

2.10.2 Set probe in electron rejection mode 

2.10.3 Observe probe output 

2.10.3.1 Probe voltage versus time 

2.10.3.2 Probe current versus voltage 

2.10.3.4 Ion temperature 

2.10.3.5 Total ion density 

2.10.4 Record all future data 
EP2.il Ion Mass Spectrometer 

2.11.1 Ion mass spectrometer to GO status 

2.11.1.1 Spectrometer to ON 

2.11.1.2 Wait for spectrometer READY 

2.11.1.3 Spectrometer to GO 

2.11.2 Set mass range 

2.11.2.1 Lower mass limit = 1 Amu (Atomic Mass 
Unit) 

2. IT. 2. 2 Upper mass limit = 40 Amu 

2.11.2.3 Resolution = 1 Amu 

2.11.3 Observe spectrometer data 

2.11.3.1 Check spectrometer angle of attack 

2.11.3.2 Display ion current versus mass 

2.11.3.3 Display ion density versus mass 

2.11.3.4 Total ion density 

2.11.4 Record all future data 


t 


I 


I 


1 


1 


7^ 


EP2.12 Neutral Mass Spectrometer 

2.12.1 Neutral mass spectrometer to 60 status 

2.12.1.1 Spectrometer to ON 

2.12.1.2 Wait for spectrometer READY 

2 . 12 .1.3 Spectrometer to GO 

2.12.2 Set mass range 

2.12.2.1 Lower mass limit = 1 Amu 

2.12.2.2 Upeer mass limit = 40 Amu 

2.12.2.3 Resolution = 1 Amu 

2.12.3 Observe spectrometer data 

2.12.3.1 Check spectrometer angle of attack 

2.12.3.2 Display neutral current versus mass 

2.12.3.3 Display neutral density versus mass 

2.12.3.4 Display total neutral density 

EP2.13 Planar Electron Trap - Retarding Retention Analyzer 

2.13.1 Planar electron trap to GO status 

2.13.1.1 Trap to ON 

2.13.1.2 Wait for trap READY 

2.13.1.3 Trap to 60 

2.13.2 Choose electron mode 

^ 2.13.3 Set energy range 

2.13.3.1 Lower energy = 1 electron volt 

2.13.3.2 Upper energy = 1000 electron volts 

2.13.3.3 Sweep linear voltage versus time 

2.13.4 Observe trap data 

2.13.4.1 Display electron current versus voltage 

2.13.4.2 Electron density versus electron energy 

2.13.5 Record all future data 
EP2.14 Planar Segmented Probe Trap 

2.14.1 Planar segmented probe to 60 status 

2.14.1.1 Probe to ON 

2.14.1.2 Wait for probe READY 

2.14.1.3 Probe to 60 

2.14.2 Probe into ion mode 

2.14.3 Observe probe data 

2.14.3.1 Probe current versus time 

2.14.3.2 Display of angle of plasma flow 

2.14.4 Record all future data 


• .-f ■■ 


EP2.15 


EP2.16 

EP2.17 

EP2.18 

EP2.19 
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Special Displays 
2.15.1 


2.15.2 


Densiti es 

versus time 

2.15.1.1 

Electrons 

2.15.1.2 

0+ 

2.15.1.3 

0p+ 

2.15.1.4 

n6+ 

2.15.1.5 

0 

2.15.1.6 

02 

2.15.1.7 

N2 

2.15.1.8 

He 

Position as function 

2.15.2.1 

Local time 

2.15.2.2 

Latitude 

2.15.2.3 

A1 titude 

2.15.2.4 

Magnetic d 

2.15.2.5 

Universal 

Angular position 

2.15.3.1 

r 

2.15.3.2 

A 

2.15.3.3 

A 

2.15.3.4 

0A 

2.15.3.5 


2.15.3.6 

R 

2.15.3.7 

n 

2.15.3.8 

X 

2.15.3.9 

'Jj 

2.15.3.10 

Y = angle 


angle between Shuttle velocity 
vector and the platform axis 

How Call Up? (By Time, By Latitude, By Longitude?) 

Stop All Data Recording 

All Instruments to OFF 

Stow Boom A 


END EXPERIMENT 




I 


I 





(I 
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3 . 2.3 Theory and Background - EP 2.0 

3 . 2 . 3 . 1 Location of the Shuttle . Assume the Shuttle is in a circular 
orbit at altitude H with some nonzero inci^ination angle = 1 in degrees. 
Assume Shuttle moves from west to east. At some time the Shuttle will be 
moving southward when it crosses the equator. We begin orbit time at this 
instant at the descending node, i.eL, the point in space where the Shuttle 
crosses the equator moving southward. 

The local time and Greenwich niean time of the first descending node 
are set by the keyboard: 


Initial local time = LT^ in hours :mi nutes : seconds (201) 

Initial universal time = 6MT^ in hours:minutes:seconds (202) 

The initial east longitude of first descending node is: 

L0N„ = I (LT^H - UT^H) '^=9 (2«3) 

where LT and GMT are equal to LT H and UT H expressed in hours, 
o o ^ o o'^ 

The initial orbit time: 

OT^ = 0 seconds (20^) 


Once the Shuttleis started in its orbit the elapsed orbit time and 
the Greenwich mean time will advance in real time second for second. If 
the orbit is advanced by OTSi seconds from the descending node the equa- 
tions, computed in seconds for updating the orbit, are: 


t (seconds) = OTSi 

(205) 

UTS = UT S 
0 

(206) 


where t is the time from the start of the orbit program. We assume that 
sidereal and synodic time are equal so that the local time of the descend- 
ing node is fixed. We also assume the sun is at the equinox position and 
stays there. 


I 


1 


1 


I 


T 


i 
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Define the orbital angular velocity which is a constant of the motion 
In our circular orbit: 

% / a 

(rVh) radians/sec (207) 


where 


R • earth radius = 6371 km 
e 

_ o 2 

g = surface gravity = 9.807 x 10 km/sec 
H = orbit altitude in kilometers, a constant 
The east longitude of the Shuttle portion at any time t is given 


by: 


LON = IlON^ + tan~V- )| . , . (208) 

I o \ cosojt /'modulo 36 O deg 

The south latitude of the Shuttle is given by: 

LAT = +cos ' ("sinoit sini) ~ +sin ' (+sin(ot sini) (209) 

The local time in the orbit is given by: 

in hours if t is in seconds and LON is in degrees. 

Time is displayed in the following fashion: 

21 :45:32 UT or LT 

where 21 represents the hour, 45 represents the minutes past that hour, 
and 32 represents the seconds past that minute. 

Time is kept to the nearest second. However, whenever a time is used 
in a formula it must be expressed in either decimal hours or seconds. The 
orbit period is given by: 

TORB IT - ^ in seconds (211) 


3. 2. 3*2 Shult le Att i tude . Consider the transformation needed to 
get from an Earth-based inertial frame into th' Shuttle XYZ frame. Consider 
the quas i - i ner t i a I frame X'Y'Z' whose center is on the Shuttle XYZ origin 
but where X' always points south Y' east and Z* always points radially up- 
ward. Let gamma, lambda, and delta be the three usual Eulerian angles. 



Let X'Y'Z' be called the Earth frame even though it is traveling around the 
Earth on the Shuttle orbit. The transformation that takes vectors in this 
Earth frame into vectors in the Shuttle frame is given by 

( cosAcosf-cos/\ sinfsinA 
-s i nAcosF-cos ‘.s i nl'cosA 
s I n ‘ s I n r 

It should be noted that the X'Y'Z' frame is not a real inertial frame. 
The Shuttle attitude control s/stem must be active in order to maintain a 
fixed setting of gamma, lambda, and delta throughout an orbit. The 
X'Y'Z' coordinates are, in effect, spinning at the rate of one revolution 
per orbit about an axis perpendicular to the orbit plane passing through 
the origin of X'Y'Z'. 


cos' s i n! +cos’ cosr" s i nA s i n' s i nA, 

-5 i nAs i nF+cos ■ cosTcos.' cosAsinA ) (212) 
-sin\cosF cos' 


I 
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3>2«3>3 Shuttle Velocity Vector . In Earth centered R, 0 , 
coordinates starting at line of nodes: 

R ■ (Rg + H) 10000 = cos ^(-slnwt slni ) ■ tan ^ (tanuit cosi) 

The Shuttle velocity has components VO and V(j) only, and because of 
our choice of X', Y', Z', one can easily see that the Shuttle velocity 
vector is given by: 


\ + Vy. y- 

A A 

where X' and Y' are unit vectors and: 


Vy , = 0)R 


sin coscot 


7, . 2 ,vl/2 

(1 - sin I sin (ot) 


(213) 


Vyl = (OR 


cos I 


(1 - sin^i sin^(ot)'^^ 


Vg = coR meters/sec 


R = (R^ + H) 1000 meters 


Given this formulation of V^, one can transform it into any relevant 
coordinate axis on the Shuttle. In particular, one parameter that appears 
often is the angle between the Shuttle velocity vector and the Xj coor- 
dinate axis of the platform: 


ANGLE VX 


cos (Vg-x^) 


(oR 


(214) 


'^S * X 1 “ ^A->P ^S->A ^E-^-S I '^Y » 

\o 


(215) 


'x^ component 



where 


cosij/cosx”CosJ2sinxsin(J> cosij^sinx+cosficosxsinij^ sinfisimjj 


■ I “sInij;cosx“CosJ2s1nxcosip 


sin!i2sinx 


”Sini|/sinx+cosJ2cosxcosij> sinQcosi^ I (216) 
-sinSlcopx cosi[2 J 




cos6^cos<j)^ 

sin6^cos(j)^ 


cose^^sint})^ 


si ri(|>^sl n9^ 


-sine 


The angles 0^, 4^^, 4, x» ^ were defined in EP 1.0. 


Expl ici tly: 


j = ANGLE VXj = cos. 


^A->P ^S-^A ^e->S 


( s i n i coscotV I 

. /| 


2 2 

(1 - sin i sin wt) 


component 


3.2. 3.4 Earth's Magnetic Field . The Earth's magnetic field is 
assumed to be a centered, due-soath pointing dipole field whose dipole 
moment = 8.07 x 10^^ Teslas - m^ = M. 

The equation for the magnetic field vector in the X', Y'* Z' system 


where 


“earth = “x' * “z' 


„ 8.07 X 10 ^ X sinX 

By, — -r— in 


(R^ + H) 


gammas 


B,, = - 


1.6l4 X 10 X cosX , 


(Rg + H) 


in gammas 


I 




In this equation, X is given by: 


X = LAT + ir/Z = cos (-sinwt sini) 


0 < X <L 180 deg 


The magnetic dip angle is given by: 


»i i — , 

.. , cos B*Z‘ “1 -2cosX 




sX 

* 


= cos 


2sinu)t sini 

(3sin^u)t sin^i + 1 ) ^ 


The dip angle is measured from the local vertical. 

The total magnetic field in gammas is given by (R^ and H in 
kilometers): 

1/2 

1C 2 2 * f ^ 

i« I _ (8.07 X. 10 ^)(3sin cot sin i + 1) in gammas (222 

The magnetic field measured by the three-axis fluxgate magnetometer 
is given by: 


BX = T^_ I 0 
e^PI 


component 


\“Z'/»X 2 component 


T I 0 
e->»pl 


871/1 

^ ^ ' X^ component 
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where and is given by Equation (219) and: 

" ^A-vp ^S-»A ^e-vs (223-5) 

3-2. 3*5 Neutral Atmosphere Model 

a) Temperature . For this model assume that the neutral, ion, and 
electron temperature are equal at all times in the given altitude range of 
120 to 400 km. Inactuality, the electron temperature is usually signifi- 
cantly higher than the other two, but this fact is neglected in the model. 

Consider the temperature at a given altitude H, latitude LAT, and 
local time LT. All temperatures will be in °K. The maximum temperature 
at the al ti tude is : 

TMAX = 2100 - 9000 (224) 

The temperature at some LAT and LT at this altitude H is given by: 

T = + 0.28 sin’*^|LAT/2|) (l + A cos^*^(J)) 

where 

A = (0-^8)[*^os^’^IlAT/ 2| - sin^‘^[LAT/2|] 

1 + 0.28 sin’ '^1 lAT/ 2| \ 

We want the sine and cosine of the absolute value of LAT/2 to the 1.5th 
power. 

T is an angle that depends upon the solar hour angle, SHA: 

t = SHA - 45 deg + 12 deg sin(SHA +45 deg) (227) 

The solar angle is given in terms of the local time expressed in decimal 
hours: 

SHA = ||80 deg * ii (360 deg)| (228) 

Thus given H, LAT, and LT, the temperature at the Shuttle's position is 
defined by Equation (225)- T will always represent this temperature. 


(225) 

(226) 


.1 


I, 


I 


b) Neutral Atmosphere Density . For this model, consider the 
neutral atmosphere to consist of only four components: atomic oxygen, 

molecular oxygen, molecular nitrogen, and helium. The model will be 
defined in terms of the maximum density of these constituents at the 
120 km altitude level: 


NMAXq = 7-60 X 10 atoms per cubic meter (2 

NMAXq 2 " 7.50 X 10^^ molecules per cubic meter 
NMAXj ^2 “ ^*00 ^ 10 '^ molecules per cubic meter' 

1 3 

NMAXj^g = 3*^0 X 10 atoms per cubic meter 
(U.S. Standard Atmosphere, 1966) 

Each constituent will be distributed as a function of altitude and tem- 
perature or altitude and LAT and LT by way of Equation (225) • 

Define a scale height for each species: 


kT(Rg + H)‘ 

M.g R ^ 
e 


M.w (R + H) 
I e 


where 


= atomic weight in AMU species i 

= temperature at the point of measurement 
in OK. 

= 0.880 km/°K 


Therefore; 


0.880 T(Rg + H)‘ 


(231) 


8i* 




«02 = 


0.880 T (Rg + H)‘ 
32 R^^ 


H 


0.880 T (R + H)‘ 


e 


«HE- 


0.880 T (R. + H)' 

c 

k 


Each species is distributed in the foi lowing fashion: 


dn. 


i _ ^ ^ 

H. 


( 232 ) 


SO that: 


N(H, LAT, LT). = NMAX. {■^)exp 


120 H, 


(233) 


gives the density of species i, at altitude H, latitude LAT, and local 
time LT. If the variation of (R^ + H) with H is neglected because 
R >> H, then the integral can be done and: 


N. = NMAX 
1 



M.R 


0.880T (R+H) (0,0137) 

CO e 


M.R (120 - H) 


exp- 


I e 


0.880T„(Rg+H)' 


( 234 ) 


where 


T<„ = (j^)(l + 0.28 sin';^|LAT/2l)(l + A cos^*^’ j) 


: J-5 


2.5 T, 
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Refer to Equations (224) and (225) for H = “. Tj 2 q is given by Equation 
(225) with H f 120 km. NMAX. is given for each species, i, in Equation 
(229). M. is the atomic weight in AMU of species i: 

Mfl - 16 

”02 “ 32 

M - 28 
"2 

M - 4 
HE 


» f j:5 


Each neutral species density is now fully described as a function of 
altitude, local time, and latitude or of altitude, orbit inclination, 
and orbit time. The total neutral density at any point in the orbit is 
given by: 


NTOTAL = N(0) + N(02) + N(N2) + N(HE) 


(235) 


3.2. 3.6 Ionospheric Density Model (Electrons) . The equation that 
governs the density of electrons or ions in the ionosphere is: 


dn 2 

q - an for electrons 


(236) 


where 

q = the photo ionization rate 
a = the recombination rate 

It has been assumed that all ions have equal recombination rates so 
that the ratio of any single ion species to the total ion density does 
not change. 

a) Electron Density . From Nawrocki and Papa, "Atmospheric 
Processes," page 1-61, an empirical formula for the photoionization rate 
is obtained. We put in a latitude dependence. If iLATj < 85 degrees: 


120 < H < 180 km q. = 3000 cosJlAT 


(237.1) 


\ ll 


J80 < H < AOO km q = 300O cos|LAT|exp 


180 “ H 


(237.2) 


If LAT > 85 degrees: 


120 < H < 180 km = 3000 cos [85 deg| 


(238.1) 


180 < H < ^00 km q^ = 3000 cosi85 degjexp ~ ~ (238.2) 

Thus, there are four separate regions where q is defined so that one must 
determine LAT and H before defining the proper q value to be used below. 

Therefore, at a fixed H and LAT, the local time dependence of the 
electron density is given by Equation (236). We have simplified consider- 
ably by making q a constant during daylight. Of course at night q = 0 
and we assume night begins at all altitudes and latitudes at I8OO LT and 
ends at O6OO LT. 

During the day O6OO 5 LT < I8OO. Equation (236) has the solution. 

LT is measured in decimal hours: 


where 


= tanh 


n(LT) - nnanh[a^'i|U^* C,] 

tanh fl X 10'* .') . , 

. -1 \ n“ / 

3nh - T v 7 Tv 

/, , q .A^ 3 2_ . , lp ..V 10 \ 


(239) 


At night the solution must be divided into the two i ntervals 1 8OO < 
LT < 2AOO and 000 < LT <' O6OO because of the way we keep LT < 2400. 

For 1800 < LT < 2400: 


n(LT) 


*fq (LT - 18) 3600 . T 

= n + 

>- n J 
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where 




-1 


For 0000 < LT < 0600: 


n(LT) 


*r<i (LT + 

= n[_- 


6) j 1 oi ^ 3600 

modulo 2k 


+ C. 


-1 


(2^1) 


n' is the daytime equilibrium density at latitude LAT and H. It is 
related to the equatorial equilibrium concentration by: 


and 


1 /2 

n “ cos (LAT) for 0 deg < |LAT| < 85 deg 


n = n^ cos'^^ 85 deg for ILAT] >85 deg 


(2li2.l) 


(242.2) 


All that remains to be done is to define the altitude dependence of 
n^. An empirical relationship for 120 <H < 240: 

= 4.9 X 10^^ exp .—-^^9) electrons/m^ (243.1) 

For 240 < H < 400: 

n^ = 4.9 X 10 '^ exp ") electrons/m^ (243.2) 

This completes the description of the electron density as a function of 
altitude, local time, and latitude. 

b) Ion Density . This model assumes that only the ions of atomic 
oxygen, molecular oxygen, and nitric oxide are present. The masses of 
these ions are : 

O'^ = 16 AMU (244) 

O^"*" = 32 AMU 
no'*' = 30 AMU 
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Assume that at any given altitude the ion density of a single species 
is directly proportional to the electron density at that altitude: 

N(02+) = F(02+) n(LT) (245) 

N(N0+) * F(N0+) n(LT) 

N(0+) = [I ’ F(02+) - F(N0+)3 n(LT) 


where F(02+) is this fraction for 02 ^, F(N0+) for NO^, and n(LT) is given 
by Equations (239) to (241). 

The fractions are defined empirically in three altitude ranges. For 
120 < H < 210 km: 


F(02+) - (0.436)exp(y^ “ (246.1) 

F(N 0 +) = (o.8l4)exp(-j^ - -j^) (246.2) 

For 210 < H < 240: 

F'(02+) = (0.22)exp(9.57 - ^ (247.1) 

F(N0+) = (0.47)exp(8.63 ~ (247-2) 

Finally, for 240 < H: 

F(02+) = (0.22)exp(6.05 - ^+ 155 ) (248.1) 

F(N0+) = (0.47)exp(5.11 " ^ + -j^) (248.2) 


Each of the three species' densities is now fully defined as a function of 
altitude, local time, and latitude. 


' ■ - ’ ■— - - i - - i I J I 
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3 . 2. 3 . 7 Instruments 

a) Electric Field Receiver and Magnetic Field Receiver : In this simu- 

lation these two instruments are turned on, but they respond only to the 
random noise source and there Is no way of calculating the mean square field 
value. A request for the plot of mean square field should be answered with 
NOT AVAILABLE FOR EP 2.0. Otherwise the two Instruments are fully functional 
using the random noise generator as a signal source. 

B) Fluxgate Magnetometer and Rubidium Magnetometer . Equations (219) 
through (223) give all needed information to determine the magnetic fields 
observed by these instruments. 

C) The Cylindrical Electron Probe . The electron temperature is given 
by Equations (224) to (228). 

The electron density is given by Equations (237) to (243). 

The current versus voltage display must be calculated in the following 
fashion: Let 

4 ,i P = probe potential in volts 

T = electron temperature in °K 

"3 

n(LT) = electron density in (meter) ^ Equations (239) to (241) 

"19 

e = electron charge = 1.6 x 10 ^ coulombs 

"23 o 

k = Boltzmann's constant = 1.38 x 10 joules/°K 

-31 

m = electron mass = 9.11 x 10 kg 

A = probe area of cyl I nder diameter 0.1 cm, 35 cm long = 

-3 2 

1.1 X 10 meter 

V = Shuttle velocity in METERS/SEC not km/sec 

s 

= Platform unit vector In X 2 direction. 

There are two cases for the probe current. If P £ 0, probe current is 
given by 

1 " = I - I. I electron current (249) 

e « e 

^ [!■ ... 

I . ion current 

I 



.1 


I 


I 






Remember that P is negative 


iT » A n(LT) e 




'^s * ^2 “ ^A-+P ^S-»-A ^e-^ ''y* 


V 

1 

V 


component 


NOTE: is in meters/sec. 

DefinecI .1 as positive when electrons flow into the probe. 
If P > 0 


= l"^ - it 

e I 


,+ „ A n(LT)_e |^' . ' | 

TT ' S 2' 
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(250) 


( 251 ) 


( 252 ) 


(253) 

( 25 ^) 


1^ = A n (LT) e 
e 


+ (eP/m) 


1/2 


[(s) 

(' ' •(/?))] 


exp (-eP/kT) 


where 


(255) 




dt = error function. 


(4 


The angle between the probe axis and the velocity vector is given by 


CPV = cos 


-1 ^s * ^2^ 
Rw 


(256) 


O) Spherical Probe . Each one of the ions 0^, NO^, 0^ contributes a 
separate current to the spherical probe. Each one of the three ion currents 
has the form given approximately by 


. 2 „ s , 

I . = irr N . e -g- 1 




i/2 -M.(V;-V)‘ 

I e s 

— m 


(257) 


r^ is the sphere's radius =0.05 meter 

N. is the density of ion species i given by Equation {2k5) 

"19 

e = 1.6 X 10 coulombs 

Vg = shuttle velocity = Rw x 10 see Equation (213)* Vg must be 
measured in meters/second 


2 r "t 

$(x) = — / e dt = error function 
0 ■ 


M. = particular ion mass in kilograms 


Mq+ = 16 X 1.66 x 10 'kg 


M«+ = 32 X 1 .66 X 10 ' kg 

^^2 


(258) 


= 30 X 1.66 X lO"^' kg 


V = (2dP/M.) 

e I ' 


( 259 ) 


P = electric potential on the probe >0 in volts 

T = temperature Equation (225) in °K 

"23 o 

k = Boltzmann constant = 1.38 x 10 ^ joules/ K 


I 


i 
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The total ion current at any voltage, P, is given by the sum of each of the 
currents at that voltage. 


* measured ^ ^ ^NO'*’ 


( 260 ) 


The ion temperature is given by 


T. = T 
ion 


(261) 


T is given by Equation (225). 

Ion density is given by 

Ion density = n(LT) is given in Equations (239), (2^0), and (2i»l) 
depending on the local time, altitude and latitude. 

( 262 ) 

The plot of log^Q 'measured simply what it says: 


log base ten of probe current nanoamps 


+9 + log^Q 'measured 


Current in nanoamps = 10^ I j of Equation (260) 

measured 


(263) 


(26ii) 


E) Planar Segmented Probe . For this experiment we use the probe to 
investigate the angle between its axis and the ion streaming velocity 


Angle PSPV = Angle VX, as given by Equation (2l8) 


(265) 


'■W 


To shorten the notation below, let 


Y = Angle PSPV 


( 266 ) 


The probe is constructed with radius 0,05 meter and the sides of the probe 
are 0.05 meter high. For this simulation it is not necessary to find the 
current In each probe segment. The total current is sufficient. The expres- 
sion for the current is in two parts. If 


I 


i 


I 


tan y < 2 


I « 2r^ eNR w cos y j^sin ' (l - (l " ^ 


tan Y " tan 


(267) 


r = 0.05 meter 
e = 1 . 6 X 1 0 ^ ^ c 


oulombs 


N =» total ion density = N(LT) given by Equations (239) to (2^1) 
R = Orb iter geocentric radius in meters 
w = orbital angular velocity in Equation (207) 

Y = angle PSPV in Equations (266) or (218) 


I f tan Y >2, then 


(268) 

(269) 


Current measured in nanoamps = 10^1 of Equation (267) (269) 

Log base 10 of probe current in nanoamps = 9 * (270) 

F) Planar Electron Trap . The current from the planar trap at any 
potential P, assuming that the ions are being rejected Is 


1/2 

|(P) =>r^(-^) ej^N(LT)/f exp(-mV^^/2kT) + /T^ exp(-mV^^/2kT^ ) 


+ N„/t 7 exp (-mV V2kT, 


( 271 ) 


where 


r = 0.05 meter 
k = r .38 X lo"^^ joules/°K 

- I 'Q 

e = 1.6 x 10 coulombs 

N(LT) = electron density given by Equations (239) to (2i^l) 
T = temperature given by Equation (225) 




m » 9.11 X 10“^^ kg 

8 7 . 

- 7.16 X 10 electrons/m'^ 
k o 

= 2 X 10 K 

N 2 = ^.75 X 10 electrons/m^ 

T., = i».6li X 10® °K 
= 2eP/m 

P *» trap potential in volts. 

Log base 10 of current in nanoamps = 9 + (272) 

Log base 10 of electron energy in electron volts = ^ (273) 


The log base 10 of the electron density in electrons per cubic meter is the 
log of the density of electrons in the range E to E = dE; E = energy. 



exp(-eP/kT) 



exp(-eP/kT^) 



exp(-eP/kT2) 


I 


(27^) 


All terms defined as in Equation (271). 

Log base 10 of electron density = Equation (27l») 

(275) 

G) Ion Mass Spectrometer . For the ion mass spectrometer there will 
be data only at masses 0 = 16, NO = 30, and O 2 = 32. 

The log base 10 of ion density in ions per meter‘d is for ion Mass = 16 
log base 10 ion density = logjQN(O^) from Equation (245) (276) 


log^QN(NO^) from Equation (245) (277) 


log^QN(02) from Equation (245) (278) 


Ion Mass = 30 AMU 

log base 10 iori density = 
Ion Mass =32 AMU 

log base 10 ion density = 


AMU 
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x>- 


All other masses have no display of log base 10 of ion density. 
A general expression for the ion current is 


where 


l(Mj) • irr^ Nj e cos y[l + ^(Vg cos y »¥j//2kT) 


+ (2kT/M|)'^^ exp-(Mj cos^y/2kT)(\/^cos 


r » 0.02 meter 

N. is the number density of species i equals Equation (2A5) 
~ 19 

e » 1.6 X 10 coulombs 


Vg * Shuttle velocity in meters/second 


cos y *» OTg • x^/lVgl Equations (211») to (218) 


2 r 

$(x) = — / e dt = error function 
■^0 

M. =: mass of each species in kg. See Equation (258) 

T = temperature See Equation (225) 

, k = 1.38 X 10 joules/°K 

The current at mass 16> 30 and 32 is obtained by choosing the proper 
M. in this above equation. All other masses have zero current. 

Q 

Ion current in nanoamps = 10^ I (M.) of Equation (279) 

Total ion density = N(LT) of Equations (239) to (2l»l) 

The angle between spectrometer axis and Shuttle velocity vector is 


ANGLE IMSV = ANGLE VX^ of Equations (2lA) or (218) 


(279) 


N. and 

I 

( 280 ) 

(281) 

( 282 ) 
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H) Neutral Mass Spectrometer . The angle between spectrometer axis and 
velocity vector is the same as Ion mass spectrometer in Equation (282). 

The total neutral density is given by Equation (235) (283) 

This model has only neutral masses k, 16, 28, 32 AMU, all other mass 
numbers have no displays for Neutral Mass - k 

log base 10 density = logjQN(HE) from Equation (23^) for HE (28A) 
Neutral Mass =16 

log base 10 density = log^QN(O) from Equation (23^) for 0 (285) 


Neutral Mass = 28 


log base 10 density = log^QN(N2) from Equation (234) for N2 (286) 


/i Is 


Neutral Mass = 32 


log base 10 density = log^QN(02) from Equation (234) for 02 (287) 

The neutral current is given by Equation (279) where for 
Mass 4 

M. = 4 X 1.6 X lo”^^ kg 

N. = N(HE) from Equation (234) for HE 


Mass 16 


M. = 16 X 1.6 X lo"^^ kg 

N. = N(0) from Equation (234) for 0 


( 289 ) 


a I- 


Mass 28 


M. = 28 X 1.6 X lo'^^ kg 

I 


N. = N(N 2 ) from Equation (234) for N2 


Mass 32 

M. = 32 X 1 .6 X lo"^^ kg 

N. = N(02) from Equation (23i|) for 02 

, I i 

all other currents are zero. 

Q 

Ion current in nanoamps = 10^ l(M.) from Equation (279) 

3.2. 3. 8 Special Displays . 

A) Ion densities are given by Equation (2A5) 

B) Electron density is given' by Equations (239) to (2Al) 

C) Neutral densities are given by Equation (23^) 

D) Temperature is given by Equat ion (225) 

E) Times are explained in Equations (201) to (210) 

F) Locations are explained in Equations (201) to (210) 

G) Magnetic field is explained in Equations (219) to (223) 

H) All the angles except XIV are set by the keyboard. Angle 
XIV is given by Equation (218). 

THE END 
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(291) 


(292) 
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3.3 IONOSPHERIC MEASUREMENTS WITH THE SUBSATELLITE 

3 . 3.1 Experiment Description 

Ionospheric Measurements with the Subsatellite EP 3.0. 

The subsatellite system which provides the ability to make measurements 
at great distances from the Shuttle is an extremely valuable tool for atmos- 
pheric and space physics measurements. In this simulation the procedures 
for launching the subsatellite are outlined in a general fashion, and a 
large number of the subsatellite instruments are used. The experiment is 
concerned with the measurement of Ionospheric densities, species and tem- 
perature at large distances from the Shuttle. The ionospheric model is the 
same one used in EP 2.0. The altitude. Inclination and attitude of the 
subsatellite can be varied over a wide range of values. 


3 . 3.2 Experiment Procedure 

EP 3*0 Ionospheric Measurements with the Subsatellite 
EP 3.1 Deploy Subsatellite 

4.1.1 Check area clear for ejection 

4. 1 . 1 . 1 TV picture 

4. 1.1. 2 Safety and warning systems 

4.1.2 Eject Subsatellite 

4. 1.2.1 Check subsatellite power systems GO 

4. 1.2. 2 Check subsatellite telemetry system GO 

4. 1 . 2.3 Check subsatellite transponder system GO 

4. 1.2. 4 Check subsatellite control system GO 

4. 1.2. 5 If all systems GO push ejection button 

4 . 1.3 Turn subsatellite instruments on 

4. 1 . 3 . 1 3“axis fluxgate magnetometer 

4. 1 . 3.2 Cylindrical electron probe 

4. 1 . 3 . 3 Segmented plana trap 

4. 1 . 3.4 Ion mass spectrometer 

4. 1 . 3.5 Electric field meter 

EP 3*2 Determine Subsatell i te Position and Attitude 

4.2.1 Set subsatellite altitude = 300 km 

4.2.2 Set subsatel 1 i te orbit inclination =10® 

4 . 2.3 Set local time of descending node = 0000:00 

4.2.4 Set universal time of descending node =0000:00 

4 . 2.5 Set spin rate at one revolution per orbit 
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V, 

HjT 


3 . ' 





^.2.6 Set initial attitude angles 

i*.2.5.1 SPHl = 90° - is = 80° 
k.2.5.2 STHETA = 270° 

it. 2.5. 3 SPSI = 0° 

k.2.7 Start subsatel 1 i te orbi t 
EP 3-3 Generate Subsatel iite Data 

4.3*1 Check 3"axls magnetometer outputs 

4 . 3.2 Check cylindrical probe outputs 

4 . 3.3 Check segmented trap outputs 

4 . 3.4 Check ion mass spectrometer output 

4 . 3.5 Check electric field meter outputs 

4 . 3.6 Record all data 

EP 3.4 Turn All Subsatellite Instruments Off 

4.4.1 3"axis fluxgate off 

4.4.2 Cylindrical probe off 

4 . 4.3 Segmented trap off 

4.4.4 Ion mass spectrometer off 

4 . 4.5 Electric field meter off 

3 . 3.3 Theory and Background: Ionospheric Measurements 

with the Subsatellite 

3.3.3. 1 Subsatellite Orbit . The subsatellite orbit is treated exactly 
1 ike the Orbi ter. 

3 . 3 . 3 . 2 Subsatellite Attitude . The three axes of the subsatellite 

are denoted by s^, s^, s^. Y^, represents the coordinate system 

fixed on the subsatel 1 i te orbit with X always pointing southward, Y 

always pointing eastward, and always pointed upward. Let the subsatellite 
velocity in the x’, y\ Z^ system be V^^IX^ + Vy1y' for an orb i t incl ined 
at i degrees at altitude H then 


V = R 


s i n i cos wt 


w ,, .2. . 2 .]/2 

(1 - sin i sin wt) 


cos I 


w ,, .2. .2 ,xl/2 

(I - sin I sin wt) 


AvA 


/ 


(312) 


100 


where R = (R + H) x 1000 and Vyl and VI are in meters/sec, and 

6 AT 


R / q \ ^ 

Where w^ = — -rj- — +^1 “ subsatellite orbital angular velocity 

V® J 

R = 6371 km 
e 

g = 9.807 X lO"^ 

H altitude in kilometers 

It is assumed that the subsatellite will rotate in a positive direc- 
tion about its own s^ axis with an angular velocity '^^p* 


Wgp = 2 tt (subsatellite spin rate) 


(313) 


If we apply the transformation 




x->v 


VI VI 

-^cos(w^p-w)t 


VI VI 

w 'WT 


“x' 


cos (w -w) t 
sp 


sin(w^^-w)t 

sp 


-sin(w -w)t 
sp 


-cos (w -w) t 
sp 



w 


to a vector in the 


(314) 


we transform it into a frame whose z-axis is 
perpendicular to the orbit plane and which is rotating about its z-axis 


with angular velocity w 


sp 


This is exactly the frame in which we wish to 


place the subsatellite axes s^ 


s^ in order to make the measurements 


of the physical phenomena 

in the direction of V if w - w 

sp 


With this orientation the s^ axis always points 
THEREFORE WE WILL ALWAYS PUT THE SUB- 


SATELLITE IN ITS ORBIT WITH s^ perpendicular to the orbit plane and s^ 
aligned with V at the beginning of the orbit (t =0). If w then 
s^ will always stay aligned. ^s->^V will then transform to; 
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/V V oA 


s-*-V 




•Vyl V^l 


( 315 ) 


y 


If not, then will rotate about out of the local horizontal. 

Define the three Eulerian angles PHI, THETA, PS1 which rotate , yV, 
z' into Sj, s^, Sy In general the transformation 


= /cos Tj) COS <() "- cos 0 sin (j> sin ij) cos sin f|) + cos 9 cos <j) s i n tj) 

“sIn cos (f> - cos 6 sin (|) cos ip -sin ip sin (f) + cos 6 cos 4> cos ip 

\s i n 0 s i n (f) -s i n 0 cos (ji 

Rotates a vector from xN^z' into s^S2S^. In particular if 


sin ij) sin 0^ 
cos Ip sin 0 I 
cos 0 


y 


<j) = cos ^ - - 2 *^ * j ' 7 ' 2 » ® and ip = (w^ - w) t (316) 

(1 - sin k sin wt) ^ 

then SjS^Sj will be the frame described I^Y Equation (^1^) " the one we 

want to make our observation from. The initial values of ((), 0, ifi will always 
be 


Initial (j) = 4 = n/2 - I 

T 


_ 


Ini f ial 0 = 0 = -tt/2 „ 

o Z 


Initial 4-4 = 0 

o 


(317) 

(318) 

(319) 
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3. 3.3.3 Display of Subsatellite Realtime Values 

3.3.3. 3.1 Orbit and Attitude . Universal time In hours is given by 
UT = UT^H + UT^H in decimal hours is converted from GMT^ (keyboard 

entry) in hours, units, seconds. 

Magnetic dip angle is given by 

.. _ _ “1 2 sin wt sin i 

dip angle = cos = -s ^ (320) 

(3 sin^ wt sin^ i + 1) 


PHI is given by 

PHI . cos-' .. sin I cos^wt (jj,) 

(1 - s 1 n'^ I 5 in^ wt) 

THETA is given by 

THETA = 270° (322) 


PSl is given by 


PSl = (w - w) t 
sp 


(323) 


The procedure for EP 3-0 is written with w^^ = w so that 

PSl = 0° 


(32A) 


The line of sight distance to Shuttle is given by Equation (311). 
3. 3 .3. 3. 2 Subsatel 1 i te FI uxgate Magnetometer . 


in gammas 


1 - 8.07 X 10^^ sin X 
" ■ (R 

e 

„ y “1 .61 A X 10^^ cos X . 

B 1 = in gammas 

^ (R^ + H)3 


(326) 


A = LAT + Tr/2 - cos \ (-sin wt sin I) 
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For subsatellite fluxgate whose three axis X, Y, Z correspond to Sj, 
s^t Sj, respectively, we have 


(328) 


(329) 


(330) 


(331) 

where is given by (4l4.5) and , B^l by Equations (325 and 326). 

3. 3. 3. 3. 3 Subsatellite Cylindrical Electron Probe . Remembering that 
everywhere a local time, altitude, orbit time or velocity is mentioned it 
refers to the subsatel 1 i te and not to the Shuttle we can take over the whole 
formalism of EP 2.0. 

Electron temperature is given by Equations (22^) to (228) EP 2.0. 

Electron density is given by Equations (237) to (243) EP 2.0. 

Current versus voltage displays is given by Equations (249)to (255) 
of EP 2.0. 

The angle between probe axis and the velocity vector Is given by 

CPV = ir/2 + (Wgp - w) t (332) 

3.3-3. 3-^ Subsatel 1 i te Planar Segmented Trap . The current is given 
by Equations (267) to (270) EP 2.0 except that the angle y is given by 

Y = (w^„ - w) t (333) 

op 


BTOT = 


(8.07 X 10^^) (3 sin^ wt sin^ i + 1)^^^ 

(K + H)^ 


in gammas 


- ■'x-sv I “ , 

Bzl 


s^ component 


BY = T 


x->sv 


\V/ 


component 


° '^x-sv 1 ® , 

h' 


Sj component 


^ ’ 


4 


[i 


i. 


10 ^ 

3*3»3*3*5 Subsatellite Ion Mass Spectrometer . Equations (276) through 
(281) describe the ion mass spectrometer current and density displays. 

The angle between spectrometer axis and the velocity vector is: 

y = ANGLE IMSV = (w^p - w) t (33^) 

3 "3.3 *3 *6 Subsatellite Electric Field Meter , for this simulation the 
only electric field will be that produced by the motion of the subsatellite 
across the magnetic field lines. 

The electric field is given by 

Southward “ ^x' ’ <®z'> mi I llvol ts/meter (335) 

^eastward " millivolts/meter (336) 

^upward “ ” (VyO (B^^l) (10 ^) mi II ivolts/meter (337) 

and B^l and Byl are given exactly by Equations (325) and (326). The total 
electric field in millivolts/meter is 

ETOT = (E^l + Eyl + E^l)^'^^ (338) 


.f-V 
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3. A ELECTRON ACCELERATOR BEAM MEASUREMENTS 

3.4.1 Experiment Description 

Electron Accelerator Beam Measurements EP 4.0, 

Electron beam experiments have already been done in the ionosphere and 
the technique is applicable in many areas of magnetospher i c and plasmas 
physics in space. The AMPS Shuttle payloads will probably contain a high 
powered electron accelerator. 

Before the electron beams can be used in ionospheric experiments from 
the Shuttle it will be necessary to map out the characteristics of the 
beam as It leaves the accelerator. The following simulation attempts to 
describe the initial behavior of the beam as it expands from its initial 
dimensions because of its own self-field and because of the Earth's magnetic 
field. 

The beam current, energy, direction of propagation and divergence 
are controllable by the experimenter over a wide range of values. The 
measurements of the beam are made by a Faraday cup placed seven meters 
from the accelerator's exit aperture. In this simulation the magnetic 
steering coils of the accelerator are used to sweep the beam back and 
forth across the face of the Faraday cup so that the beam profile can be 
mapped without moving the Faraday cup. 


I 


I 


I 
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3.k.Z Experiment Procedure 

EP ^.0 Electron Accelerator Beam Measurement 

EP ^».l Determine Shuttle Orbit 

^1.1.1 Set Shuttle altitude at AOO km 
A. 1.2 Set Shuttle orbit inclination at 0° 

, ■ i*.1.3 Set local time of descending node = 0000:00 

* A. 1.4 Set universal time of descending node = 0000:00 

EP 4.2 Set Shuttle attitude with Z-axis parallel to Earth's 
Magnetic field 

4.2.1 Set = 270° 

4.2.2 Set A = 90° 

4.2.3 Set A = 0° 

4.2.4 Hold Shuttle in this attitude for remainder 

of experiment 

EP 4.3 Set up accelerator lens focus 

4.3.1 Set diverting lens voltage = 100 volts 

4.3.2 Set converging lens voltage = 100 volts 

EP 4.4 Set up accelerator scanning 

1 ! 4.4.1 Set BPHl =0° 

4.4.2 Set cathode heater current = 4.0 amps 

4.4.3 Set beam into programmed scanning mode 

4 . 4 . 3 . 1 Selected programmed nx5de 

4. 4 . 3 . 2 Set Wj = 2 tt 

4. 4. 3 . 3 Set W 2 - 20ir 

EP 4.5 Accelerator tests 

4 . 5.1 Test beam azimuth and deflection action 

4 . 5.2 Test beam current value 

4 . 5.3 Test grad current value 

EP 4.6 Turn Faraday cup on 

’ 4.6.1 Set photoelectron suppressor voltage 

4.6.2 Set Faraday cup current amplifier gain 

4 . 6.3 Set Faraday cup bias and scale factor 

EP 4.7 Final . accelerator preparation 

4 . 7.1 Set accelerator voltage P to 10,000 volts 

4 . 7.2 Set cathode heater current to 4,0 amps 

4 . 7.3 Set grid voltage to 100 volts 

4 . 7.4 Set firing duration to 1.00 seconds 

4 . 7.5 Turn charging current ON, observe accelerator 
READY 

EP 4.8 Time-to-go until firing 

• ■ ' r ' ■ ' ■ 

4.8.1 Ser time-’to-go at 60 seconds 
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A0 


I 


(I t: 


EP 4.9 Set countdown to GO wait until beam has fired 
EP 4.10 Check data 

4.10.1 Check electron accelerator beam current 
versus time 

4.10.2 Check beam firing angles 

4 . 10.3 Check Faraday cup current 

4.10.4 Check beam angle with respect to Earth’s 
magnetic field 

EP 4.11 Ref i re 

4.11.1 If want to refire beam 

4.11.1.1 No changes in beam parameters 

4.11.1.1 reset firing time in 4.7 
and continue 

4.11.1.2 Change beam parameters 4.11.1.2 
Repeat steps 4,4 to 4.10 with 
desired changes 

4 . 11.2 If no ref i re 

4.11.2.1 Turn accelerator OFF 

4.11.2.2 Turn Faraday cup OFF 

3 . 4.3 Theory of Background 

3 . 4 . 3 . 1 Motion of an Electron in an Electron Beam Fired in the 

Presence of a Magnetic Field (General Information). In the 

following an attempt has been made to describe the current that would be 

measured by a Faraday cup placed a fixed distance above the exit aperture 

of an electron accelerator. The Faraday cup has an entrance aperture area 

of 10 square centimeters. The exit aperture of the electron accelerator 

is circular of radius 50 centimeters. 

An attempt is made to describe the Faraday cup current as a function 
of: 

1) Accelerator current (l) 

2) Accelerator Potential = Beam Energy (P) 

3) Orientation of the ejected beam with respect to the ambient 
magnetic field (fipcjij) 

4) The magnitude of the ambient magnetic field (B) 

5) The distance between the Faraday cup and the accelerator 
aperture (Z^) 

6) The beam dispersion angle (ifi) 
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The electron accelerator is located on the Shuttle X-axis, a distance 
from the origin of the Shuttle X, Y, Z coordinate system. It is assumed 
that the Faraday cup is mounted at a distance = 700 centimeters from the 
accelerator. 

3*^*3.2 The Beam Radius. The general motion of the particle in the 
beam is given by 


m dv _ ^ V X B \ 

y-E = 4 ttp 


V-B 
V X E 


= 0 


V X B = 


li ^ 

C 6t 


hjL j + 1 

C C 6t 


<Sp 

6t 


+ V*J = 0 


(^ 01 ) 


subject to the initial position and velocity of the electron and the beam. 

7 1/2 

The initial density p^ = 1.3^ x 10 IP I = accelerator current in 
amps. P = accelerator potential in volts and p^ is in electrons per cubic 
centimeter. 


[A closed form, self-consistent solution for the electron position 
cannot be generated and so it is necessary to make many approximations in 
order to bring the problem into relative simplicity. We shall divide the 
problem; into two parts. The first treats the expansion of the beam due 
to the self-electric field. The second treats the expansion due to the 
beam divergence. 

3 . 4 . 3 . 2.1 Radial Expansion due to Self-Field . For this problem let 
us go into a cylindrical coordinate system (r, <}>, Z) that travels up the 
Z-axis at the beam velocity. The equation of motion are: 




109 


** *2 eE r • 

r - = — jjp + rw| 


r = 


iL 

dt 


r(t> + 2f(j) = wr (j) * 


_ d4» 


dt 


(^ 02 ) 

(403) 


w ^ 


2 2 

Q _ 2e p ttR 
eB eE o o 


me m 


mr 


(404) 


= radius of beam when expansion starts = 50 cm 


“28 

M = electron mass =9-1 x 10 grams 


v-10 


e = electric charge of electron = 4.8 x 10 ESU 

7 “ 1/2 

= electron density at t = 0 = 1.34 x 10 IP el/cc 
c “ speed of light = 3 x 10^^ m/sec 

We have at t = 0, <|) = 0, <^ = 0, r = R^, f = 0 as the initial conditions 


Equation (403) gives ^ = y\“T” " ^ 

' r 


(405) 


and Equation (402) gives 


2 4 

2 2 w^R ^ 

a w r ■ o 


(< 106 ) 


integrating 


£ 

/T 


— = la in r/R 



(407) 


We shall assume in this section 





(408) 
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so that 


R 

K _ O 


I ^ I 


/uTiVR 


dx 
a/ Imx 




o 2 

u j 

e du 


where r is now the radius to which the beam has expanded after a time 
The integral cannot be done in closed form so vye expand it roughly by 


2 / 

f Vi ^ / 1 + e^^ 

I e dx = I — ;; ib or 


« (1 + r/R) /In r/R 


= (1 + r/R) -QO (r/R^ - l) 


o' 7 


in the region 1 r/R^ _< 10 approximately. This gives approximately 


or 


(-' * //C)= sf 


t/Rq = (1 + 

''o 


R = 50 cm 
o 


a = 1.46 X 10® R cm/sec 

o 


Thus the pure electric field expansion will be approximated by 


rjj = 50 (1 +10.22 X 10® cm 


(409) 

(410) 


(411) 


(412) 


(413) 


(414) 


(415) 


in 


where is the radius of the beam at time t. 

All of the many approximations given above will prove to be not too 
bad so long as the measurements are made within 10 meters of the accelerator 
aperture. 

3 . ^.3 .2. 2 Radial Expansion due to Beam Dispersion 


D I . riT. 

R; - — sin (♦ sinj^-' 

b 

7 1/2 

= beam velocity = 5-93 x 10' P cm/sec 
cB 7 *“ 1 

w = — = 1.76 X lO' B sec B = mag field in gauss 
me 

= dispersion angle 


(416) 


r, = beam radius at time t 

b 

r = distance from center of beam < r. 

— b 

After a time t an electron which started at distance r will travel 

o 

to a distance 


rj (t) = [1 + 2x^ (1 - X 1/x^ cos (wt + ctn ^ x))]^^^ r^ (417) 
R' 

where x = — 
r 

o 

If the density of the beam were uniformly distributed along the 
r-direction from 0 to r^^ then the electron charge density at time t would 
be given by 

p R ^ 

p(t) =-^5-2- (418) 

However because the actual dispersion of the beam is given by 

ip sin j— (419) 

b 

instead of ij; alone, the density is not uniformly distributed in r. 
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The actual density becomes 


p(r,t) = 


o 2 2 (ti> sin ■r“) 

R ‘ IT cos Zri 

o o b 


r/(t) 4/ - / (1 + ttV4) 


(420) 


provided is not too large, if •" £ r^^ and 


p(r ,t) H 0 for r > r. 


(421) 


3. 4. 3. 2. 3 Combined Dispersion and Self-Field Radius 

We shall assume that the electric field acts upon a dispersively 
expanding beam as follows: 


r^ (t) = r^ [1 + 2 X (1 - /j ^ j^^2(cos wt + ctn x)] 


1 aJ/2 


(422) 


where r' is the self field radius 
o 


r = r [1 + 10.22 X 10^ t]^^^ 


o o 


(423) 


In particular the total radius of the beam after a time t is given 


by 


r. (t) = 50 [1 + 1,02 X 10^ 1^^^ [1 + 2 x^ (1 - /l + l/x^ 


”1 1 /2 

cos (wt + Ctn x) ] cm 




1 = beam current in amps 

P - beam potential in volts 
"2 1/2 ”1 

x = 6.75 X 10 (P B sin i|^), P in volts, B in gauss 

w = — =1.76 X 10^ B sec ' 
me 
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and p(r,t) is given by Equation (420) and (421) with «'j^(t) given by (424 

and R = 50 cm. 
o 

3. 4. 3. 3 Motion of the Center of the Beam. If the Shuttle Z-axis is 
not pointed along the magnetic field direction then 

Z 



the angles 6 and y define its position and the transformation wi 1 1 

take a vector in the Shuttle X, Y, Z system into a vector in the magnetic 
field system Xj , , Z^ which has the Z^ axis pointed in the direction 

of the magnetic field 

I I A 

e = |B| z, 


whe re 


B 




-sin Y cos Y 0 

-cos 6 cos Y -cos 5 sin Y sin 6 

sin 6 cos Y sin 6 sin y cos 5 


and 


sin Y = 


cos 6 = 


X y 


B 


B ^ + B ^ + B ^)’^^ 
z X y 


(425) 

The beam is ejected with velocity V|^ at angles BOMEGA and BPHI in the 
Shuttle X, Y, Z system 

Let, ini tial ly 


I .# % 

BOMEGA = n 

(i*26) 

BPHI = <f> 

o 


In the magnetic system Xj , Zj 



let the orientation angels be fij and c|)^ where 

cos n. = cos Q cos 6 - sin sin 6 cos (6 -y) (^27) 

loo o 

and 

(f)^ = ctn ^ [sin 6 ctn csc((|)^ -y) " cos 6 ctn ((j)^ -y)] (428) 

The motion of the beam in the magnetic field is helical; the sum of a 

uniform motion in the Z, direction at velocity V, cos and a circular 

I b I 


motion in the X,-Y, plane with velocity V, sin 12, at frequency w = — . 
Looking at the projection of the motion of the beam in the X^^-Y^ plane 
we see 



The center of the circle Is at X = R sin d),, Y, = cos (f>,. 

I c I I c I 

Let us go into a coordinate system 


X2 “ Xj + R^ sin (}.j 


V. sin n, 
D I 


Yj = Y, - R_. cos 

which is centered at the center of the circle, then the center of the beam 
moves with coordinates 

- R^ sin(wt + <j>p 


y- = R_ cos(wt + <j),) 



H6 





I I 1 


and finally let us consider the coordinates system of the beam center 

Xj = sin (wt + <|)^) 

(431) 

We see that if 

+ y3^ 1 Tj^^(t) (432) 

then we are in the beam and will measure an electron flux. If 

^3^ + V3^ > (^33) 

we are outside of the beam and will not measure a flux. Going back from 
the Xj, y^ to the Shuttle X, Y, Z system we have 

2 2 
r = [sin Y (Xjtj - X) + Y cos Y + (sin 4)^ - sin (wt + <j)^))] + 

[cos 6 cos Y (X|^ “ X) “ Y cos 6sinY'*'Zsin6 + R^ (434) 

2 

(cos(wt + <j>|) " cos 4>^)J 

where Xj^ is the X-coordinate of the electron accelerator. We arbitrarily 
take 

Xj^ = 200 cm (435) 

This value of r’ is the one for which p(r,t) is defined in Equations 
(420) and (421). 

Finally we shall need to know the velocity of the beam in theZ- 
di recti on 


dt 


[sin 6 sin sin (wt + 4>j) + cos 6 cos 


(436) 


3.i».3.^ Electron Accelerator Storage Bank . The storage bank supplies 
the energy for the electron beam. Ordinarily when an experiment is started 
the energy stored will be zero and the experimenter will ask the computer 
to charge the storage bank. We as&ume the charging current comes from a 
500 volt source and that the maximum energy stored is 100 kilojoules. 

Let 

E(t) be energy stored at time t in Joules 

i . be charging current in amps 
ch 

D(t ) = time to reach full charge 
b 

E(0) = energy stored at time zero assumed = 0 joules 

E - maximum energy stored = 100,000 joules 
o 

D(tjj) = 17.3 In (100 - 0.32 \/eT^ (^37) 

If E(0) = 0 then * 80. This represents an approximate value 

of the number of seconds to reach full charge. The charging current is 

'ch " (^3 “ ^^2?: ' 7 ) ®^P(-0-058 t) (^38) 

from a 500 volt source. 

There are two cases for E(tj^), one when the beam has been fired and 
some energy remains and the other with E(tj^) starting from zero energy. 

^.^.1 Beam Energy Remaining from Previous Firing 

Let T|^ be the time duration that the beam was fired on the 
previous firing. Let be the energy stored at the time the beam started 
firing the previous time, then the energy left after the previous firing, 
taking into account that the charging current was left on during the firing 
is 

E(tjj) = -IP^ + (0.4) [500 + (v/z.SE^ - 500) exp -0.058 Tjj]^ ( 439 ) 

I - beam current in amps P= beam potential in volts 
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3.4.3.^.! Beam Not Being Fired 

If the beam has never been fired and the bank is charging up from 
0 joules. 

E(t) = 100,000 (1 - exp - 0.058t)^ (41,0) 

E(tj^) = E^(l - exp - 0.058t)^ (44l) 

actually E(t) ^ E^ for finite tg but if we round E(tj^) upward to E^ at 
E(tg) =* (E^ - 2000) joules we obtain realistic charging times. 

If the beam has been fired for a time x^, first calculate E(t,) from 

o b 

Eq. ^39. Then calculate E(tg), from: 


E(tg) = 0.4 


500 + is/l.S E(tj^) - 500) exp - 0.058 tg 


(442) 


3. 4. 3. 5 Accelerator Current . The accelerator current depends upon 
the temperature of the cathode and the potential of the grid. The 
temperature of the cathode depends upon the cathode heater current. 

Vg = grid voltage 

1^ " cathode heater current '^h ” cathode heater voltage 

I g = beam current 

1^ and Vg are set by the keyboard. 


I = grid current 
9 
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Ig * ( 7.05 + 75 Vg exp(-n.A Vg l^■’''^)amps (A43) 

and is independent of t as long as beam is being fired. 


1 B 1 0 ^ 1 B m i 1 1 i amps 
9 0 


(444) 

V. = 61, volts 
h h 


(445) 

^.I».3.6 Storage Bank Drain Time. 

The beam power is 

I^P and we calcu 

late the approximate drain time of 

the bank by 


E 



Xj=Dra,n Time- p . , 

D 

seconds 

(446) 

1 given by Eq. 443 




If drain time >1000 seconds or If IgP £ 11500 watts then the display 
should read 

Drain Time = >1000 seconds. 

3 • ^ • 3 . 7 Cathode Heater . 

Cathode heater current = set from keyboard. 

Heater voltage t 6lj^ 

3 . 4 . 3 . 8 Diverging and Converging Lens . The lens voltages are set 
from the keyboard and control the divergence of the beam by 

m) 

= converging lens voltage. 


Tj; = ( 0 . 09 ) exp 


v^-v 




where V, = diverging lens voltage and V 
d c 


3.^.3 '9 Beam Aiming Coils . The angles BOMEGA = and BPHI = 
are set from the keyboard. If the BPHI and BOMEGA scan program is chosen 
the computer will compute values of and so that 

= w, im 

a 

= -^ (1 + sin t^) (4A9) 

where t^ starts at t = 0 when the beam is fired and stops when t = Tj^. The 
display values of Beam Firing Angles are taken from Eqs. ^^8 and 

The coll currents are 

2 cos (J). 

X-axis coil current = —f/o (^50) 

P tan 

P = accelerator potential in volts 

2 sin 

Y-axis coil current = — (^51) 

P’''^ tan 

D 

3.^.3 *10 Grid Voltage and Current . Grid voltage is set from the 
keyboard. Equation ^i»3 gives the test beam current (identical to beam 
current) and Equation gives the grid current. 

3 . ^ . 3 • 1 1 Accelerator Potential is Set from the Keyboard . 

3.4.3.12 Firing the Beam . It is important that the beam is not 
allowed to fire longer than the approximate drain time given by Eq. 446. 

3.4.3.13 Beam Angle with Respect to the Earth's Magnetic Field . I f 
the Shuttle orbit and attitude are allowed to vary then in the Shuttle 
X,Y,Z system 








By = T ^c 
X e->S 




\®z- 


X component 
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(^♦52) 




By/ 


^ = Vs ' ° 


B, = T - 
Z e->S 


®zv 

l\‘\ 

0 

\h'l 


Y component 


Z component 


(453) 


(454) 


where , B^, are given in gammas by Eq. 219 and is given by Eq. 212. 

Equation 425 gives the expression for y and 6 in terms of the B^^ By B^ 
of Eqs. 452-54. 


The angle between the beam axis and the magnetic field is given by 
of Eq. 427 except we must now allow and (|>^ to vary if the beam scan 
program Is used, i.e., replace by and c|)^ by 


cos = cos cos 5 - sin sin 6 cos((j)j^ - y) 


(456) 


where SJ^^and <})|^ are given by Eqs. 448 and 449. 

3.4.3.14 The Faraday Cup Current . The electron current into Faraday 
cup is 

I = e (area of cup) (density of beam at cup) (ve loci ty normal to cup) 
■ c 


= (10)(e)(p(r,t)) 


dz 

dt 
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First we must determine if the cup is in the beam. The time between 
an electron leaving the accelerator and its reaching the cup is 


t^ = ^ (sin 6 sin sin <|»^ + cos 6 cos ' seconds 


(^♦ 58 ) 


where Z 700 cm is, the cup d'istance above the accelerator and V. is the 
c I/O ** 

beam velocity = 5-93 x 10' x P cm/ sec, and and (()j given by Eqs. it27 and 428. 

1) There will be no current to the cup for timer t < t^ 

2) After time tj the current to the cup will be if r r^ 


I = 1.27 X 10 
cup 


^ |cos^iJ< s i n j i - i|;^(l + 


-1 


(459) 


[sin 6 sin sin(wt^ + (j>j) + cos 6 cos n^jx 


-1 


[1 + x^(l ’ Vl + l/x^ cos (wt, + Ctn"^x)l X 


where 


[1 + 1.02 X 10^ t^l amps. 


I, = accelerator current Eq. 443 
b 


P = accelerator potential 
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■%Jr 


= 2R ^(1-cos wt,) + Z ^ sin^ 5 + 2Z R sin 6 [cos (wt, + (j),) - cos <|).] 

C iC CC li i 


Vl sin 1/0 -I 4 .C 

R = i- = 3.37 P |Bt‘ I sin SI. 10 ^ cm 

c w ' I ’ I 


(460) 


Z = 700 cm 
c 


2500 [1 + 2x^ (1 - + l/x^)cos (wt^ + ctn 'x)] 

[1 + 1.02 X 10 ^ 

b 1 

X = 6.75 X lO"^ iB^^l sin ip x 10^ 


(461) 


(461.1) 


ip from Equation (447) 

6 from Equation (425) 
from Equation (456) 

<|)^ from Equation (428) with and (j>^ 

w = — - = 1.76 X 10^ \bJ X lO"^ secs' 
t^ from Equation (458) 


({)|^ of (448) and (449) 


-1 


If r > r. then 
b 


I = 0 amps 
cup 


(cup in amps) ~ * (cup in mi 11 lamps) 


(462) 


The only time dependence of the cup after t = t^ comes from the fact 
that and <()^ will vary in time of the beam scan program is used. 


Jt .-V 
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3.5 LIDAR TRACE OF ACOUSTIC GRAVITY WAVES IN THE SODIUM LAYER 

3. 5 - 1 Experiment Descriptions EP 5-0 Lidar Trace of Acoustical Gravity 

Waves in the Sodium Layer 

Description 

Atmospheric sodium is concentrated in a layer from 70 to 100 km above 
the ground. The characteristics of this sodium layer have been studied both 
experimentally and theoretically. Recent advances in laser technology have 
made possible accurate measurements of sodium densities versus altitude. 

The AMPS Shuttle payload can provide a suitable platform for these lidar 
measurements of the sodium layer. 

One experiment that is possible, if the return signals from the sodium 
layer are large enough, is the determination of the existence of acoustic 
gravity waves traveling in and through the sodium layer. In the following 
simulation the magnitude of the density variations caused by the gravity 
waves is modeled from experimental data and then the ^ImMl^ted laser data 
is used to reproduce the modeled wave densities. The detection of the waves 
has been simplified by assuming that all of the wave energy is in one fre- 
quency mode and that the unperturbed sodium density is constant throughout 
the 80-100 kilometer altitude range. ■ i 

The experimenter is able to control the laser pulse energy, wavelength 
and the pulse repetition rate. Preanalysis of the return laser signals into 
specified altitude intervals has been assumed so that the actual dat^ rates 
are not beyond the capability of the CVT hardware. During ail actual AMPS 
flight, this kind of preanalysis will be done by a special lidar dedicated 
computer. 

3 . 5 - 2 Ej^periment Procedures EP 5-0 Lidar Trace of Acoustical Cifavity 

Waves in the Sodium Layer 

EP 5 -i Determine Shuttle Orbit 

5 - 1.1 Set Shuttle altitude at I80 km 

5.1.2 Set Shuttle orbit inclination at 27 ° 

5.1.3 Set local time of descending mode = 1800:00 

5.1. A Set universal time of descending mode » 1800:00 

EP 5.2 Set Shuttle al t i tude wi th z-axi s toward nadir 

5.2.1 Set T = 180° 

5.2.2 Set A =180° 








'-'Jjis' 
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5.2.3 Set A = 0 ° 

5. 2 . A Hold Shuttle in this position for remainder of experiment 

EP 5.3 Set up Ildar system 

5.3.1 Select dye 

5.3.1 .1 Rhodamine 6 G 

5.3.2 Select spectral control elements 

5.3.2 . 1 Broad tune grating 

5. 3. 2. 2 Fine tuning Fabry-Perot 

5.3.3 Select laser output monitors 

5.3.3 . 1 Power meter 

5. 3. 3 . 2 Absorption cell 

5. 3.3 . 3 Spectrometer 

5. 3. 3. A Photometer #3 “ Direct 

5. 3. A Select power supply control 

5. 3 . ^.1 Manual fire 

5 . 3 . ^.2 Select maximum power = 1000 watts 

5. 3*^. 3 Select laser pulse energy 0.1 joules 

5.3. ^.^ Select laser pulse repetition rate = 2 Hz 

5.3.5 Adjust laser transmitter 

5. 3. 5.1 Select automatic alignment 

5.3.5 . 2 Select beam divergence = 1.00 milliradians 
5. 3. 5’. 3 Select azimuth = 0 . 00 ° 

5.3.5.^ Select elevation = 0 . 00 ° 

5 . 3 . <i Adjust 1 idar receiver 1 

5. 3 * 6*1 Set aperture =100 

5. 3. 6.2 Set convergence = 0 

5. 3 . 6. 3 Select sodium D interference filter 

5. 3.6.4 Select Fabry-Perot filter 

5. 3. 6 . 5 Set detector high voltage at 3000 volts 

5. 3.6.6 Select detector 

5. 3. 6 . 7 Set field stop = 1 .00 mi 11 iradians 

EP 5.4 Activate 1 idar system 

5.4.1 Open dust cover 

5.4.2 Turn lidar dye pump to GO 

5.4.2 . 1 Check that dye I s flowing 
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5.4.3 Put broad tune grating in the beam 

5.4.4 Select broad tune wavelength = 589.0 nm 

5.4.5 Turn Fabry-Perot tuning motor to ON 

5.4.6 Set fine tune to 0 

5.4.7 Put Fabry-Perot in the beam 

5.4.8 Set wavelength resolution at 0.01 nm 

5.4.9 Turn spectrometer on lidar output to 60 
5.4.9.! Turn photometer #\ to GO 

5. 4. 9. 2 Set upper wavelength range = 8OO.O nm 

5. 4. 9-3 Set lower wavelength range = 300.0 nm 

5. 4. 9. 4 Set resolution = 10 nm 

5.4.10 Turn power meter to GO 

5.4.10.1 Set upper wavelength = 590 nm 

5.4.10.2 Set lower wavelength = 58O nm 

5.4.11 Turn phototube #3"direct to 60 
5*4.12 Adjust absorption cell 

5*4.12.1 Put absorption cell in the beam 

5.4.12.2 Turn cell temperature control to 60 

5.4.12.3 Set cell temperature- 300°C 

5.4.12.4 Turn phototube #2 to 60 
5.4.13 Turn lidar flashlamp power supply to GO 
5*4.14 Set maximum power level 1000 watts 
5.4.15 Select fire time as fire at wi 11 

EP 5* 5 Fire laser 

,5.5.1 Turn laser fire to GO (laser will fire at 2 Hz, 0.1 joules 
until laser fire is turned to NO-GO) 


EP 5*6 Adjust laser spectral output 

5.6.1 Observe display of spectrometer on laser output 

5.6.2 Adjust fine tuning Fabry-Perot until line width is 0.01 nm 
and central wavelength is 588.996 nm 

5.6.3 Check absorption cell transmission =0% 

EP 5*7 Stop laser select laser fire = NO-GO 

EP 5*8 Select laser energy per pulse = 10 joules 

EP 5.9 Fire laser select laser fire = GO 






EP 5.10 Check laser output 

5 . 10.1 Read output power = 10 joules 

5 .10 .2 Read absorption cell transmission = 0% 

EP 5.11 Stop laser fire; laser fire = NO-GO 
EP 5.12 Select automatic laser fire 

EP 5.13 Activate transmitter 

5 . 13.1 Automatic alignment system to GO 

5 . 13.2 Beam divergence to GO 

5.13- 2.1 Check beam divergence = l.OO mrad 
5-13-3 Azimuth and elevation to GO 

5-13-3-1 Check azimuth = 0.00° 

5.13- 3-2 Check elevation = 0 . 00 ° 

EP 5-l‘^ Activate receiver 

5- 1^-1 Turn detector high voltage power supply to GO 

5.1 A. 2 Turn aperture control to GO 
5 .IA .3 Turn field stop to GO 

5.14.3.1 Check field stop = 100.00 mrad 

5.14.4 Put sodium D filter in the beam 
5-14.5 Adjust Fabry-Perot filter 

5 . 14 . 5.1 Adjustment motor to ON 

5-14.-5-2 Set central wavelength = 588.99 nm 
5-14.5-3 Set bandpass = +0.01 nm 

EP 5. 15 Select data analysis parameters 

5 . 15.1 Set upper altitude = 110 km 
5-15-2 Set lower altitude =70 km 
5-15-3 Set altitude resolution = 0.10 km 

EP 5.16 Fire laser; select laser fire = GO 

EP 5-17 Check Ildar data summary display 
5-17-1 Shuttle altitude = 180 km 

5.17- 2 Shuttle latitude 

5.17- 3 Shuttle longitude 
5-17-4 Shuttle local time 

5 . 17.5 Sun is down = YES (If sun down = NO, discontinue experiment 
until sun is down = YES) 
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5. 17.6 Angle between z-axis and nadir = 0° 

5.17.7 Laser power = 10 joules 

5.17.8 Laser central wavelength = 588.99 nm 

5.17.9 Laser divergence = 1 . 00 mi 1 1 i rad ians 

5.17.10 Receiver field stop = 100.00 milliradians 

5.17.11 Receiver central wavelength = 588.99 nm 

5.17.12 Data upper altitude = 110 km 

5.17.13 Data lower altitude = 70 km 
5.17.1^ Data height resolution = 0,10 km 

EP 5.18 Check data displays 

5.18.1 Raw count rate versus altitude 

5.18.1.1 If count rate less than 1000 counts adjust laser 
parameters to higher counts 

5.18.2 Sodium density versus altitude 

5.18.3 Average sodium density versus altitude 

5. 18. it Sodium density minus average density versus altitude 

5.18.5 Sodium density minus average sodium density versus position 

5.18.5.1 Select altitude = 90 km 

5.18.5.2 Select altitude = 100 km 

5.18.5.3 Select altitude = 80 km 

5.18.6 Sodium density minus average sodium density versus altitude 
versus position (3“dimensional projected display) 

EP 5.19 Gather data for one event in local time 

5.19.1 Record all data and instrument settings 

EP 5*20 Stop laser fire; select laser fire = NO-GO 

EP 5*21 Turn instruments off 

5.21.1 Dust cover to closed 

5.21.2 Fabry- Perot filter adjustment motor to OFF 

5.21.3 Sodium D filter out of beam 
5. 21. A Field stop to OFF 

5.21.5 Aperture control to OFF 

5.21.6 Detector high voltage to OFF 

5.21.7 Azimuth and elevation to OFF 

5.21.8 Beam divergence to OFF 

5.21.9 Automatic alignment system to OFF 


5 . 21.10 Lidar flashlamp power supply to OFF 

5 . 21.11 Phototube #2 to OFF 

5 . 21.12 Set absorption cell temperature to = 30 °C 

5.21.13 Absorption cell out of beam 
5. 21.1 A Phototube # 3 *’direct to OFF 

5.21.15 Power meter to OFF 

5.21 .16 Phototube §\ to OFF 

5.21.17 Spectrometer on lidar output to OFF 

5.21.18 Fine tuning Fabry-Perot out of beam 

5.21.19 Fabry-Perot tuning motor to OFF 

5.21.20 Broad band tune grating out of beam 

5.21.21 Absorption cell temperature control to OFF 

5.21.22 Lidar dye pump to OFF 


END OF EXPERIMENT 
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3.5.3 Theory and Background: EP 5.0 LIdar Trace of Acoustical Gravity 

Waves in the Sodium Layer 

3.5.3. 1 Wave Theory 

According to Yeh and Liu [Rev. Geophys. and Sp. Sc. J_2, 193 (197^)1 
the dispersion relation for the acoustic gravity wave assuming an isothermal 
atmosphere is given by: 


(k ^ + k 

X y 




+ k 2 = 

2 

c 

o 


('■'■“7) 


(501) 


and the group velocity is given by 


^ r- , / 2 2\ ^ / 2 2. - , 21/ 2,, k 2,, 2^. 2v 2\ 

Vg = j^x k^(o) -tOj^ ) + y ky(w ) + z k^cu / (u (k^ +k^ )c^ j 


if 


From Yeh and Liu we can obtain values for 


2 % „ , irt "2 -1 

“b “ 300 ' 


c “ 280 m/sec = speed of sound 
o 


(502) 

(503) 

(504) 


in the altitude range from 80 to 100 kilometers. 

For our purposes here we can take k^ = 0. Kochanski [J.G.R. 69., 3651 
(1964)] has shown that at the 8O-IOO km range the major wave-like components 
of the atmosphere are traveling mainly to the East with a velocity of about 
40 meters/sec and that the Vertical wavelength is 10 kilometers. We shall 
call the x~direction East. Then if 


“fc 

k = a3 — 
x c 


0) = Ot 0). 


(505) 


the dispersion relation gives 
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1.23 - ° o 




3" = 


( 506 ) 


1 - a‘ 


and the velocity expression gives 


1 - 


V 3 

X 


(507) 


1 - 


^ 

3c„ 


Using the experimental data of Kochanski , one obtains 


= 5.84 X 10 ^ meters ^ 


V = 40 meters/sec 

X 


( 508 ) 

(509) 


Using these values in Eqs. 506 and 507 one can obtain the approximate 
solut ions 


3 = 10.5 


(510) 


and 


a = 0.5 


(5M) 


These give 


-h <*1 

k = 3.94 X 10 meters 


(512) 


0<jT 

X = 7 -—= 16 kilometers 
X k . 


U) = 


1.05 X 10 ^ sec ^ 


T = 600 sec = period of wave, 
w 


(514) 

(515) 
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For our purposes we shall consider that all of the acoustic gravity 
waves are restricted to this one mode in the altitude range 80-100 km. We 
shall try to observe these waves by measuring the sodium density as a func- 
tion of z « altitude and x - eastward position. 

To obtain the density fluctuation associated the wave we shall 

equate the time-averaged energy density with the wave's kinetic energy den- 
sity. Yeh and Liu give 


E = 


1/2 






^ 22,2 2v 

2oi c (w - 0 ). 1 
o b 


( 516 ) 


where p* is the first order term of the pressure variation. This transforms 
Into 




(a^-1) (a^-e^) 


X 2 „2 

+ Y H 



(517) 


where 


Y = l.A 

H = scale height = 5 km 
p' = first order density variation 

= zero order ambient density ot, 3» c^» k^ as before 

We set 



P V 
o X 



and obtain 


p' = 0.0319 pQ 


(518) 

(519) 

( 520 ) 
(521) 



(523) 
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Therefore, the density at any point x, z, y and time t is given by 

p(x,y,z,t) » p^(x,y,z,t)^l + 0.0319 exp i x - z - wt)j (524) 

where we have assumed the wave is traveling upward. Or since 
P ^ p(z) 

p(x,z,t) = p^j^l + 0.0319 cos(k^ X - k^ z - wt)j (525) 

where k^ is given by Eq. 512, k^ by Eq. 508 and w by Eq. 514. 

3.S.3.2 The Ambient Sodium Density . For the purposes of this simu- 
lation we shall define the sodium density as follows: 


p(Na) = 3 x 10 ^ atoms/cm^ for 80 _< z £ 100 km 

and 

p(Na) = 0 for z < 80 km 
z > 100 km 

where z is the altitude from the ground. 

3. 5. 3*3 The Wave Perturbed Sodium Density . The following function 
describes the sodium density 

80 £ z < 100 

p(Na) ” 3 X lO^fl + .0319 cos(.394x - .548(z-90) - *01 t)*j atoms/cc 

^ ( 528 ) 

where X is the distance Eastward in kilometers on a given latitude line from 
0°L0N. z is the altitude in kilometers and t is the local time in seconds 

p(Na) = 0 for 80 > z > 100 (529) 

If LON = Shuttle longitude in degrees and LAT = Shuttle latitude and 
H = Shuttle altitude in kilometers, then 


( 526 ) 


(527) 


V 


1 


I 


X = (R^ + H) (cos LAT)(L0 N)j|^ kilometers 


(530) 


R = 3371 kilometers 
e 


3. 5. 3.1* Day and Night . If 

(R + H)^ cos LAT sin [LON + 15(UT-18)] > R (H^ + (531) 

then the Shuttle is In the dark and the sodium measurement can be made. 

If the Shuttle is In sunlight It is unlikely that the return signal from 
the 80-100 km sodium layer could be detected. Also if the moon is up, It 
may make the measurement impossible depending on the spectral distribution 
of moonlight and the Earth's albedo. However, we shall Ignore lunar effects 
here. 

3.5. 3.5 Laser Pulse . The number of photons in the laser pulse Is 
given by 


Ej 6 X 10'® Xjj 
"ph ■ 1250 


(532) 


E = pulse energy In joules 


X. = laser pulse wavelength in nanometers 


3.5«3.6 Laser Interaction with Sodium Layer . Let 


Z = altitude of observed sodium above the ground 
p(Na) = sodium density (Eq. 525) 

i - ■ ■ 

e = detector efficiency =0.2 

a = sodium atomic absorption cross section 
o 

N . = number of photons in the laser pulse 
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Az = height Interval being observed 

2 

A = receiver collecting area = .785 m 
c 

then the number of photons counted by the detector from the height interval 
Az is 

N = A^e Az rt N . p(Na)[ATr(H-z)]"^ (533) 

c c o pn 

where we have assumed that the optical depth of the layer is very small. 

For a laser tuned to the sodium D line 

= 8.97 X 10 cm^ (T is sodium temperature in °K) (53^) 

X = 588.996 nm 

and the counts at this wavelength are 

. p(Na) Az E 
N =3.32x10 

p(Na) in atoms/cc 
E = laser energy in joules 
Az, H, Z in kilometers 
if 588.996 jjD.l nm then = 0 

3.5. ^ Details of Experiment Displays 

3.5. ^.! Laser Dye . Rhodamine 6G is by far the best dye to use for the 
sodium experiment. It will give a detectable response in the wavelength 
range 565“625 nm. We shall approximate its response by 

, . . (X - 588.9)^ 

relative response = exp - (537) 

i.e. if you are not at X = 588.9 nm you get less light energy. ’ 


(535) 


(536) 
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3. 5- ^.2 Dye Pump . The dye pump must be GO before: 

Laser can fire. "Dye In flowing" is automatically YES 
if pump is GO 

Number of firings with this dye is simply the running 
total of fires during a given simulation. 

3, 5. 4, 3 Laser Wavelength Settings . The laser wavelength is roughly 
set by broad tune grating and finely by the Fabry-Perot. Actual laser 
wavelength 

= BT + 0.1 exp(-^^ (538) 

where BT is the broad-tune wavelength set by keyboard and 
FT is the fine tune setting set by keyboard. 

The actual line has a shape given by 

l(X) = I exp -[(X-Xj^/RES^l (539) 

O A/ 

where RES is the wavelength resolution set by keyboard and where 1 is the 
intensity. 

3.5.^*^ Laser Output Spectra . The spectrometer on the lidar output 
measures the spectral content of the laser beam. The upper and lower 
wavelengths determine the range of the display. The resolution determines 
how many points are plotted. 

If the range is 300 to 800 nm and the resolution is 10 nm, there will 
be 50 points plotted, if resolution =0.001 nm then 5 x 10 points would be 
plotted. If the ranges were 588.8 to 589*2 and resolution = 0.001 then 
kOO points would be plotted. 
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The formula for the displayed curve of laser output spectrum Is 


I (X) = exp 


- 588.9)' 


(X - X^)' 


900 


L 


(RES)' 


(5^»0) 


where X = wavelength of plotted point 

X^ is given by Eq. 538 in nanometers. 

3, 5. <4. 5 Laser Output Energy and Photon Count . If the spectral 
range on the power meter does not Include X„ then energy should 
= 0 j ou 1 es . 

Let E^ = pulse energy set from keyboard the actual pulse energy will be 






oules 


(5^1) 


5 . 7 . 1 ^ - 

Photon content is given by 

"ph “ \ 5 X photons 


(5't2) 


given by Eq. 5^1 
X given by Eq. 538. 


3. 5. 4. 6 Absorption Cell . The absorption cell heats up linearly 
and stops at temperature set by keyboard 


T = (10 C)t + T 


(5^3) 


where T^ is cold cell temperature =30 C t is time in seconds from turning 
on the cell., Example: keyboard set for T = 500 C, Turn on cell takes 

t = 47 seconds to get to 500°C and then it stays fixed. 





Percent transmission Is given by 

% = 100[1 - exp ' (\^ - 588.996)^1 ( 5 ^^) 

3 . 5*^*7 Laser Firing . Power level must be set at least 100 times 
greater than the keyboard set energy per laser pulse. 

"Fire at will" means that once laser fire Is 60 the laser will fire at 
whatever rep rate Is set. If time is picked instead of "Fire at will", the 
laser will start firing at that time If all systems are GO and continue at 
the set rep rate. 

Laser flashlamps are READY if Lidar flashlamp power supply (5~17) 
is GO, 

3 . 5.^. 8 Lidar Receiver Setup . Automatic alignment must be GO to 
obtain receiver signals. Divergence is given by keyboard setting azimuth 
and elevation must be GO and must be set to 0.00 to obtain s receiver signal 
in this simulation. 

Detector high voltage must be GO and >2000 volts to have a received 
signal. Voltage level Is set by keyboard. Aperture control must be GO. 
Receiver signal depends on setting 


(— V 

\ioo; 


(5^5) 


where Nj. I s received counts and AP is aperture setting. 

Field stop must be GO, and if beam divergence is greater than the 
field stop 


< I 

C \bdI 


f FS < BD 


(5A6) 
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FS = field stop value 
BD = beam divergence 
Sodium filter must be in the beam. 

Fabry-Perot adjustment motor must be GO. The receiver response 

« exp[-(X^- Xpp)^/RESpp] (5^7) 

Xpp is the receiver Fabry-Perot w: ^length setting. •'ESpp = receiver bandpass 

3. 5-^ .9 Lidar Receiver Data Dispiays . For this simulation there is 
no detector response for altitudes out of the range 80 to 100 kilometers. 

The altitude resolution is picked to make the detector count rate 
acceptably large. The larger the resolution the more counts. 

Shuttle altitude is set by keyboard LAT given by Eq. 209; LON given by 
EQ. 208; and LT given by Eq. 210 in hours. Angle between Shuttle Z axis and 
nadir = A- 1 80*^ with A set by keyboard. All other values are set by keyboard. 

The number of counts versus altitude is a function of 

sodium density given by Eqs. 528, 529 where 
t = (LT in hours) (3600) seconds, 
altitude being observed in km 
altitude of Shuttle in km 
distance Eastward given by Eq. 530 
laser wavelength in nanometers given by Eq, 538 

data height resolution in km 
wavelength resolution of the laser beam 
wavelength resolution of the receiver 

wavelength of the receiver 

receiver aperture setting 


p(Na) = 


Z = 
H = 

X = 
= 


a 

Az 

RES 

RESpp 

Xpp 

AP 





' 1 . 







l^to 

BD = beam divergence 
FS = field stop 

E = actual laser energy given by Eq. 5^1 

oCl 

The number of counts is given by 

= 3.32 X 10^ p(Na) Az E^^^ (H-Z)"^ 

=< " (“^) ^ ■ 588.996)'] (568) 

2 

the term (X - 588.996) comes from the fact that the sodium D line filter is 

X/ 

In the receiver beam. The function Q(BD,FS) Is 
BD £ FS 
if 

Q = 1 (5^9) 

BD > FS 

1=(w)' <550) 

Because of the nature of p(Na), in this simulation N s 0 if Z Is not 

c 

in the range 80 £ Z £ 100 km. For complete accuracy one should worry about 
the fact that the term Az might include some altitudes outside of the range 
80-100 km, but no corrections will be made for this effect. 

The display will show altitude from chosen lower limit to upper limit 
and there will be one data point for each Az in that interval. 
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The display of sodium density versus altitude is a plot of Eqs. ^28 
and 529. 

p(Na) = 3 X 10^[1 + 0.0319 cos (.39Ax - .589(z-90) - .Olt)] (528) 
with a data point for each Az. 

The average sodium density is averaged over time, distance x and 
altitude z. 

p (Na) = 3 X 10^ atoms/cm^ 80 £ z < 100 (551) 

p (Na) =0 80 > z > 100 (552) 

5.7.28 - 

The density of sodium minus average density is as a function of z. 

p(Na) - ;= 95.7 cos [.39^x, - .589(2-90) - .Olt,] (553) 

where x, and t, are position and time of the previous laser firing. Get on 
point for each Az. 

The density of sodium minus average density is as a function of x 

p(Na) - p (Na) = 95.7 cos[.39^x - .589(z,-90) - .Olt] (55^) 

where z, is the altitude selected and one gets an x,t point for each laser 
fire. The operator should be able to control the total number of points, 
i.e. the position scale value. 

Finally there is 

p(Na) - p(Na) = 95.7 cos[-394x - . 589(2-90) - .Olt] (555) 

as a' function of both (x,t) and z. One set of Az points for each laser fire 
(x,t) . 


3.6 DETERMINATION OF THE WAKE OF A TEST BODY 


3.6.1 Experiment Description 

As a body moves through the ionosphere it produces a wake by sweeping 
particles out of its way. Several theories and some experiments have been 
carried out in order to investigate the wakes of satellites traveling 
through the ionosphere. The AMPS Shuttle payload offers an excellent plat- 
form from which to carry out detailed wake measurements. 

In this simulation an inflatable balloon is deployed on one of the 
50 -meter booms and the plasma diagnostic package is mounted on the other 
50 -meter boom. Because of complexity of the ionospheric wake problem the 
simulated wake ignores the effect of electric and magnetic fields. 

The experimenter is able to vary the Shuttle orbit altitude, inclina- 
tion and attitude; and the relative position of the two booms. By varying 
the boom position, he is able to trace out the wake region and measure Its 
properties. 
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3.6.2 Experiment Procedure ; EP6.0 Determination of the Wake of a 
Test Body 

EP6.1 Determine Shuttle orbit 

6.1.1 Select manual control 

6.1.2 Set Shuttle altitude = 300 km 

6.1.3 Set Shuttle orbit inclination = 10° 

6.1.4 Set local time of descending node = 0000:00 

6.1.5 Set universal time of descending node = 0000:00 UT 
EP6,2 Determine Shuttle attitude 

6.2.1 Select manual control 

6.2.2 Set r = 0° 

6.2.3 Set A = 0° 

6.2.4 Set A = 0° 

6.2.5 Hold Shuttle at these angles 

EP6.3 Deploy subsatellite 

6,3' 1 Check clear area for ejection 

6.3. 1 . 1 TV coverage 

6.3. 1.2 Safety and warning systems 

6.3.2 Eject subsatellite 

6.3.2. 1 Check subsatellite power systems GO 

6, 3*2. 2 Check subsatellite telemetry system GO 

6. 3. 2. 3 Check subsatellite transponder system GO 

6. 3. 2. 4 Check subsatellite control system GO 

6. 3. 2. 5 If all systems are GO push ejection button 

6.3.3 Turn subsatellite instruments ON 

6 . 3 . 3 . 1 Three-axis fluxgate 

6. 3 . 3 . 2 Cylindrical electron probe 

6. 3 . 3 . 3 Segmented planar trap 

6 . 3 . 3.4 Ion mass spectrometer 

6. 3 . 3.5 Electric field meter 

6 . 3.4 Determine subsatel 1 1 te pos i tion and attitude 

6 . 3 . 4.1 Set subsatel lite altitude 300 km 

6 . 3 . 4.2 Set subsatel 1 ite Orb i t inclination = 10° 

6. 3 . 4. 3 Set local time of descending node = 0000:00 

6 . 3 . 4 . 4 Set universal time uf descending node = 0000:01 

6 . 3 . 4. 5 Set spin rate at one revolution per orbit 

6 . 3 . 4. 6 Set initial attitude angles 


EP6.^ 

EP6.5 


EP6.6 


EP6.7 


EP6.8 


EP6.9 


EP6. 10 


I 1 I j 


6. 3. 4. 6.1 (j) = tt/2 - i = 80° 

s 

6. 3. A. 6. 2 e = -tt/2 = 270° 

6 . 3 . ^ . 6 . 3 ijj = 0 

Start both Shuttle and subsatellite orbits at same orbit time 
Set Boom B 

6.5.1 Set length = 25 meters 

6.5.2 Set 6g = 0° 

6.5.3 Set ({)g = 0° 

Inflate target balloons 

6.6.1 Check volume clear to inflate with Boom A alignment TV 

6.6.2 Check pressure in fill bottle 

6.6.3 Fill balloon 

6 . 6 . 3 . 1 Start fill operation 

6. 6. 3.2 Monitor balloon pressure 

6.6. 3.3 Stop fill operation when pressure = 10 Newton/meter 
Generate background data from subsatellite 

6 . 7.1 Turn instruments ON 

6 . 7 . 1.1 Cylindrical electron probe 

6. 7 . 1.2 Ion mass spectrometer 

6. 7 . 1.3 Three-axis fluxgate magnetometer 
6 . 7 . 1 .^ Segmented electron trap 

6 . 7 . 1.5 Electric field meter 

6 . 7.2 Examine all data output and record 
Set Boom A 

6.8.1 Set 6^ = tan ^ .167 “ 

6.8.2 Set 4)^ = 0° 

6 . 8.3 Set length = 20 meters 
Set platform 

6.9.1 Set X =0° 

6 . 9.2 Set = 6° = tan ' .167 

6 . 9.3 Set 4 ; = 0° 

Turn Boom A diagnostic instruments ON and look at data outputs 

6.10.1 Three-axis fluxgate 

6.10.2 Cylindrical probe 

6 . 10.3 Segmented planar trap 

6.10.4 Ion mass spectrometer 
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6 . 10.5 Rubidium magnetometer 

6.10.6 Spherical probe 

6 . 10.7 Neutral mass spectrometer 

6. 10.8 Planar trap 

EP6.11 Scan diagnostic package around balloon 

6.11.1 Record output of all Boom A diagnostics continuously 

6.11.2 Extend Boom A length from 20 to 50 meters 

6 . 11.3 Increase 8^ by 5*^ 

6.11.^ Retract Boom A length from 50 to 20 meters 

6 . 11.5 Increase 0^ by another 5° 

6.11.6 Repeat 6.11.2 to 6.11.5 until 8^ ^ 60° 

6 . 11.7 Set length = 50 meters 

6.11.8 Set 8^ = cos ' 0.60 

6 . 11.9 Scan boom around balloon 

6.11.9*1 Let runs from 0° to 360° 

6 . 11 . 9*2 As runs set platform angle x = at all times 
6.11.9*3 When (j>^ = 360° stop 

6.11.10 Decrease 0^ by 10° 

6.11.11 Set boom length at 30/cos 0^ meters 

6.11.12 Repeat Steps 6,11.9 to 6.11,11 until 0^ ± 10° and stop 
EP6. 12 Stop recording data 

EP6,13 Stow Boom A , 

6 . 13*1 All Boom A instruments OFF 
6 * 13*2 Platform angles to 0° 

6 . 13.3 Length to 2 meters 

6.13.4 0 ft - ♦a “ 

EP6.14 Stow Boom B 

6, 1^.1 Deflate balloon 
6.14.2 Boom B length to 2 meters 
EP6 . 15 Subsatellite to standby 

6 . 15*1 All subsatel 1 ite diagnostic instruments OFF 
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3 >6. 3 Theory and Background 

3*6. 3>1 The Wake > This experiment is similar to EP2.0 in many ways. 
The physical phenomena that are being measured are: l) electron density, 

2) electron temperature, 3) ion species, 4) ion temperatures, 5) ion 
densities, 6) neutral species, 7) neutral temperatures, B) neutral densi- 
ties, and 9) the DC electric field. The difference betv^een this experiment 
at EP2.0 which measured the ambient plasma is that here we deploy a target 
balloon which forms a wake region behind it as the Shuttle moves through 
the ionosphere. 

A full description of the electron and ion distribution functions 
involves the solution of the coupled set of equations 
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and the boundry conditions at the balloons surface and at infinity must be 
met. In general the solution is difficult and must be done numerically. 

In the following we assume all electrons. Ions and neutrals striking the 
balloon are absorbed, and that the flow back into the depleted region 
behind the balloon is one dimensional and is controlled by the ion and 
neutral temperature of the ambient. 







The coordinate system in which the physics is most easily worked is 
cylindrical with the z axis on the line through the center of the balloon 
which contains the Shuttle velocity vector. 

The density of any ion or neutral species with mass m. and temperature 
T which has an ambient density N(m.) given by EP 2.0 is for r R = balloon 
radius and z < 0 


2 m. V, 


p(r,z',m.) = N(m.) exp - 2kt" 


the electron density is 


(r,z' ,m ) = 2 N(m.) 


m. Ve 
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m. 

ions 
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and is independent of m^. 

For r > R and/or z' > 0 all densities are ambient densities as given 
in EP 2.0. This expression for the wake Is moderately simple in this 
coordinate system, however we must transform it into the coordinate system 
in which the measurements are being made, namely the Boom A platform 
coordinates. 

Let us define a new coordinate system at the center of the balloon, 
x'y'z' such that the z' axis lies along the Shuttle velocity vector. Then 

r = + >/x‘^ + y'^ (^) 

V/e must express the Shuttle velocity vector in this system. Start 
with the velocity in the Boom B system 


Ve = V e +V e + V e 
S XX y y z z 


where 
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See Eqs. 212, 213 of EP2.0. The transformation 
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where 
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S X y 2 


applied to gives 
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i.e,, T_, ,, creates a new coordinate system whose z-axis is pointed along 
the Shuttle velocity vector. This is precisely the x' , yV, z' coordinate 
system in which Eqs. 2, 3 and ^4 are defined. If 9^, L^, 6g, ij)g, Lg 

V- are given we wish to be able to determine if the instruments on the 
platform at the end of Boom A are within the wake region defined for 
Eq. 2; i.e., is + y^^ < R and zj^ < 0 where x^, yj^, are the coor- 

dinates of the Boom A tip in the x' , y' , z' balloon wake systems and 



( cos 0g cos <|)g cos 0g sin (})g -sin 6g\ 

-sin cos 4.J 0 

sin 0g cos (j)g sin <})g cos 9g cos 9g / 

see Eqs. 11 and 12 of EP 1.0. 

and the term (Lq + 5 meters) appears because the center of the balloon is 

D 

taken 5 meters beyond the tip of the Boom B, 

If v'x'^ + y'^ ^ R and z' < 0 then everywhere the ion, neutral or 
A A A 

electron density N(m.) or appears in an equation defining the output 

of an instrument the expressions in Eqs. 2 and 3 must be used. 

! 

3. 6. 3. 2 Boom A Instruments . ; The output of every instrument used in 
EP3.0 is the same as it was in EP2.0 so long as it is not in the wake 
region. When it is in the wake region /x^ + yj^ £ < 0, then the 

ambient density goes into the wake density ef Eqs. 2 and 3 of EP3.0. 
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Example:. I = 10°, H = 300 kfh, t = 1355 seconds, 4>0 = 9g = 0 

LTq = UTq = 0000:00 LON^ = 0 LON(t) = 90° LAT(t) = i 

LT = 0622:35 A = r = A = 0 <j)^ = 270° 0^ = 21° 

LB = 25 meters V = V =0 V = (R + H)w 

a y 6 


= 36.33 meters 



•'■A ®A 

La cos - Lj - 5. 


-1 0 0 

0 0 1 

0 I 0 



11.52 


\- 11.52 

( 12 ) 


so that zi < 0 and r = 0 is _< 5 meters we are inside the wake region. Say 
that we want to know the current to the segmented planar trap given by 
Eq. 267 in EP2.0. If the platform Enlerian angles are x “ 270° ^ 

= 270 ° then 


“COS 0, 


-sin 0, 


-sin 0, 


cos 0, 


-cos 0, 


-sin 0, 


-sin 0, 


cos 0, 


angle between velocity vector and x. platform axis is Eqs. 218, 265 and 266 


of EP 2.0. 
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Y = cos T 


A->P ^SA 




and finally we can obtain from Eq. 267 of EP2.0 
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(17) 


1 = Trr^ e N (R + H)w 


( 18 ) 


At this point we note that the N appearing in Eq. 318 EP3.0 must be 
modified because we are in the wake region. Ordinarily one would use the 
expression in Eq, 239 of EP2.0. 


„(UT) tanh [H-iiSSlHlil . c,] 


(19) 


however we must now use the total ion density from Eq. 302 of EP3.0. 


N = N(0‘‘') exp(-KM ) + N(02+) exp(-KM ^ j + N(N0+) exp - (20) 

0 \ O 2 / 

where N(0"^), N(02+) and N(N0+) are given by Eqs. 2k5 through 248 of EP2.0 
and M , M are the atomic weights in kilograms and 

o'" O2 


^ " (1K52) 


2 M. 


52 / 2kT 


Vg = Shuttle velocity in meters/sec 
k = 1.38 X 10 jou1es/°K 


T = temperature given in Eq. 225 EP2.0 

The outputs of all the instruments on Boom A can be treated in a 
simi lar fashion. 
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3 .(i, 3.3 Target Balloon . The pressure-bottle pressure always starts at 

( 63 ) 


6 2 
P = 2 X 10 Newtons /me ter 
o 


Here we must keep track of the time, T, that the f i 11 gas valve is in 
the OPEN position. We assume the pressure bottle is a sphere of radius 15 cm. 


Pressure-Bottle Pressure = 5^ + 2 x 10 exp[-(.2T)] 


(6M 


Actual Pressure in the Balloon = 5^[1 " exp(-.2T)] 


( 65 ) 


while the fill valve is open. 

Pressures both read last value with fill valve closed. Deflate value 
open causes gas to leave balloon. Let t = time since deflate valve open. 
Pressure In bottle stays at last value. Pressure in balloon is 


P, , , = P, ^ , (exp -t) 

balloon last value ^ 


( 66 ) 


Example: Experimenter opens fill valve for 10 seconds. At end of this 

time experimenter closes fill valve. 


+ 2 X 1 O*" [eXp (-2)] 

pressure bottle 


Pu 11 = 5^[1 - exp (-2)] 

ba 1 1 oon k v / j 


Then he opens deflate valve for 3 seconds and closes it. 


P L = 5^ + 2 X 10° [exp (-2)] no change 

pressure bottle l r v 


( 67 ) 


'’balloon ' [1 - exp(-2)] (exp -3) 

= 5A(e"^ - e“^) 


(68) 


END OF EXPERIMENT 


